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This dissertation explores two main topics: 1) Color Transparency and quasi- 
elastic knockout reactions involving pions and p mesons; and 2) determination of the 
J/^-nucleon scattering amplitude and scattering length via J/ip electroproduction on 
the deuteron. It is shown that at the energies available at the COMPASS experiment 
at CERN, Color Transparency should be detectable in the reaction tt + A — > it + 
p + (A — 1)* (proton knockout). It is also shown that Color Transparency should be 
detectable in the electroproduction reaction 7* + A — > p + p + (A — 1)* at small Q 2 
(where Q 2 is the virtuality of the photon) but large t (4-momentum transfer squared 
to the knocked out proton), which represents an as-yet unexplored kinematic region 
in the search for CT effects in electroproduction of vector mesons. Calculations are 
also presented for the reaction 7* + D — > J ftp + p + n at JLab energies in order to 
determine the feasibility of measuring the elastic J/^-nucleon scattering amplitude 
and/or scattering length. It is found that it may be possible to measure the J/ip- 
nucleon scattering amplitude at lower energies than previous measurements, but the 
scattering length cannot be measured. 
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Chapter 1 
INTRODUCTION 

This thesis explores two main topics: 1) Color Transparency and quasi-elastic 
knockout reactions involving pions and p mesons; and 2) determination of the J/ip- 
nucleon scattering amplitude and scattering length via J/tp electroproduction on the 
deuteron. 

The thesis is organized as follows: 

Ch. 2 discusses Color Transparency (CT) and calculation of the cross-section 
and transparency for the reaction n + A — > n + p + (A — 1)* (i.e. pion scattering 
from a nucleus of nucleon number A in which a proton is knocked-out). In short, 
Color Transparency (CT) is the vanishing of Final-State Interactions (e.g. scattering 
of the knocked-out proton by other nucleons in the nucleus) in large-momentum- 
transfer elastic or quasi-elastic nuclear reactions, and is a prediction of Quantum 
Chromodynamics (QCD). Many experiments have been done in order to search for 
the predicted effects of CT; references to these experiments are provided in Ch. 2 and 
Ch. 3. The quantity that measures the influence of Final-State Interactions is called 
the "nuclear transparency" T, which is defined as the ratio of two cross-sections: 

T = ° (1.0.1) 

&PWIA 

where a is the actual measured cross-section (either total or differential) for the re- 
action occuring in a nucleus, and <7 PWI a is the corresponding cross-section calculated 
in the Plane Wave Impulse Approximation. Complete vanishing of Final-State Inter- 
actions would give T = 1. Calculations are done here both neglecting any effects of 
possible Color Transparency, and including the effects of Color Transparency. The 
calculations are performed within the Glauber model of high-energy scattering from 
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a composite target, suitably modified to account for Color Transparency effects. It is 
shown that at the energies available at the COMPASS experiment at CERN, Color 
Transparency should be detectable in this reaction. 

Ch. 3 discusses Color Transparency and vector meson electroproduction, specifi- 
cally the reaction j*+A — > p+p+(A— I)*. Electroproduction of the p provides another 
means of detecting the effects of Color Transparency. In contrast to the purely elastic 
pion scattering considered in Ch. 2, for electroproduction there are more parameters 
that may be varied, namely the virtual photon energy v and virtuality Q 2 in addition 
to t (4-momentum-transfer-squared to the knocked out proton). These quantities, as 
well as a combination of them called the coherence length, L = „, 2iy 2 , can all affect 
the observed transparency. The coherence length plays an especially important role, 
since by varying its value the nuclear transparency T will vary even in the absence of 
any Color Transparency effects. Thus to observe an actual CT effect, one must keep 
the coherence length fixed. The calculations presented here show that CT effects may 
be observed in the reaction 7* + A — > p + p + (A — 1)*. So far experiments have 
concentrated on the large-Q 2 regime, but the results presented here show that CT 
may be observable in the small-Q 2 regime, as long as t is large enough. Calculations 
are presented here neglecting CT effects and including CT effects, and are performed 
within the Glauber model of high-energy scattering from a composite target modified 
to account for particle production and to indued Color Transparency effects. 

Ch. 4 discusses the feasibilty of measuring the J/^-nucleon scattering length 
and/or scattering amplitude in a proposed experiment at JLab, in the reaction 7* + 
D — > J/ip + p + n (where D stands for deuteron). With the mass of the J/ip being 
3.097 GeV, the threshold photon energy for photoproduction on a single nucleon is 
8.2 GeV, and is thus accessible with a 12 GeV electron beam. Most of the existing 
data on J/i/j photo- and electroproduction is at much higher energy. The upcoming 
12 GeV upgrade at JLab provides the opportunity to measure J/ip production near 
threshold (31]. The motivation for the work in the first part of this chapter (Sees. 4.2 
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- 4.4) was a proposal at JLab |30j to measure the J/^- nucleon scattering length by 
the reaction 7* + D — > J ftp +p + n, where the J/ ip is produced on one nucleon in the 
deuteron and then re-scatters from the other nucleon. The reason the J/i/'-nucleon 
scattering length is of interest is that several authors have argued that a nuclear bound 
state of the J/ip may exist [321 E3]. They propose that the force between a J/ip and 
a nucleon is purely gluonic in nature, and therefore is the analogue in QCD of the 
van der Waals force in electrodynamics, since the hadrons are color neutral objects 
(analogous to electrically neutral atoms in electrodynamics). A J/ip- nucleus bound 
state would represent a state of matter different from "ordinary" nuclei, i.e. nuclei 
composed of protons and neutrons interacting by exchange of mesons, and would 
allow investigation of an aspect of QCD (namely, gluon exchange) within the nuclear 
environment different from the usual meson exchange aspects. A diagrammatic model 
of 7*+-D — > J/ip+p+n is used for the calculations presented here, and it is determined 
that the kinematic conditions in the proposed JLab experiment and the small size 
of the contribution to the cross-section from J/^-nucleon rescattering do not allow 
the scattering length to be determined. However, it may be possible to measure the 
J/^-nucleon scattering amplitude, at higher energy, in the same experiment; this is 



discussed in Sec. 4.6. The energy of the J/^-nucleon elastic rescattering is high 
enough that many partial waves will be involved, in this case, and hence it is not 
sensitive to the value of the scattering length. The energy of the J/^-nucleon elastic 
rescattering would be in a range for which J/^-nucleon elastic scattering has not been 
measured previously (it would be significantly smaller than in the only measurements 
so far performed). 
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Chapter 2 

COLOR TRANSPARENCY AND THE REACTION 

7T + A -> 7T + P + (A ~ 1)* 

2.1 Introduction 

In this chapter, the ideas of Color Transparency are introduced, and the Glauber 
theory of high-energy scattering from nuclei is used to calculate properties of the 
reaction tt+A — > -K J r p+{A— 1)* for two cases. The first case ignores any possible effects 
of Color Transparency, while the second case includes these effects. The quantity of 
most interest here is called the nuclear transparency, which is defined as the ratio of 
two cross-sections: 

T=^- (2.1.1) 

&PWIA 

where a is the actual measured cross-section for the reaction occuring in a nucleus, 
and (JpwiA is the cross-section calculated in the Plane Wave Impulse Approximation. 
In the PWIA, all interaction of incoming and outgoing particles with nucleons in the 
nucleus are neglected, except for the interaction which is responsible for the reaction 
in the first place (e.g. ir + p — >■ n + p for the reaction n + A — >■ n + p + (A — 1)*). The 
actual cross-section a includes interactions of the incoming and outgoing particles. 
These interactions will lead in general to a value for a which is smaller than vpwiA, 
and therefore T < 1. (In the above expression for T, the cross-sections can in general 
be total cross-sections or differential cross-sections.) But a remarkable prediction of 
perturbative QCD (pQCD) is that under certain kinematic condtions, the outgoing 
particles from a reaction inside a nucleus will undergo no interaction at all with 
the other nucleons, and so the nucleus will appear "transparent" to these outgoing 
particles. The requirement for this to occur is that the reaction be a very-large- 
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momentum-transfer elastic or quasi-elastic reaction. For reasons discussed in the 
following, the reaction 71 + A — > tt -\- p + (A — 1)* should be a good candidate for 
observing the effects of Color Transparency. 



This chapter is organized as follows. In Sec. 2.2 the basic ideas of Color Trans- 
parency (CT) are discussed, along with discussion of the experimental searches for 



CT which have been performed. In Sec. 2J3 the Glauber model of high-energy 
hadron-nucleus scattering is presented. The basic results of the Glauber theory 



which will be used are presented in this section. In Sec. |2.4| and |2.5| , the Glauber 
model is used to calculate the cross-section and the transparency T for the reaction 
tt + A — > 7i + p + (A — 1)* at large pion incident momentum (200 GeV, which is the 



momentum available at the COMPASS experiment at CERN). Sec. 2.6 presents our 
conclusion, which is that at the energies available at COMPASS the effects of Color 
Transparency should be very evident. 



2.2 Color Transparency 

2.2.1 Color Transparency basics 

Color Transparency is a prediction of perturbative Quantum Chromodynamics which 
asserts that when a hadron undergoes a high-momentum-transfer elastic or quasi- 
elastic reaction inside a nucleus, the outgoing hadron experiences reduced interactions 
with the nucleons of the nucleus, compared to their interaction in free-space pp. In the 
limit of very large momentum transfer, the outgoing hadron experiences no interaction 
at all with the rest of the nucleus (it passes through without interacting with any of 
the nucleons); this is termed the vanishing of Final State Interactions. Thus the 
reason for using the word "transparency": the nucleus appears transparent to the 
outgoing hadron. For example, in the quasielastic scattering of an electron from a 
nucleus accompanied by proton knockout, A(e, e'p)(A — 1), perturbative QCD predicts 
that if the momentum transfer from the electron to the proton is large enough, the 
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knocked-out proton will experience reduced interactions with the rest of the nucleons 
on its way out. For very large momentum transfer, the fast moving proton would not 
interact with the other nucleons at all. This is in contradistinction to what would 
happen if we just sent a fast moving proton impinging on a nucleus: it certainly would 
not just pass through completely unaffected. In fact in this case, if the nucleus is large 
enough (i.e. A large enough) the proton would almost certainly undergo an inelastic 
collision with a nucleon. This can be seen from the classical result that the mean- 
free-path of a particle passing through a system of scatterers is I — ^ where a is the 
total cross-section of interaction of the particle with an individual scatterer and p is 
the number density of the scatterers. For a typical nuclear density p ~ 0.2 fmS — 3), 
and proton-nucleon total cross-section a = 40 mb (for proton momentum greater 
than a few GeV), the mean-free path is / ~ 1.25 fm. Thus for a nucleus of radius 
3 fm the incident proton would have a large probability of interacting. But for the 
fast-moving proton knocked-out of a nucleus by a hard collision, pQCD predicts that 
the probability of its interacting with the other nucleons on its way out is much less, 
and zero in the limit of very large momentum transfer. The reason why this is so is 
that pQCD predicts that the outgoing "proton" is not in fact a usual proton at all, 
but instead a system of quarks in what is called a "small-size configuration", where 
the 3 valence quarks in the proton are much closer together than they usually are in 
the proton. And pQCD predicts that the cross-section of interaction of a small-size 
color singlet with another hadron decreases the closer together the quarks are. In the 
limit of zero separation (a "point-like configuration"), the cross-section of interaction 
is zero, and to such an object the nucleus appears "transparent". This is analogous 
to what happens to a classical physical electric dipole (which of course has zero net 
electric charge): the closer together the two charges are, the smaller is the force 
exerted on the dipole by any external electric fields. In the limit as the separation 
goes to zero, the net force on the dipole goes to zero. 

It's important to note that the occurrence of the small-size configuration doesn't 
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have anything to do with the proton being inside a nucleus. For the free-space elastic 
reaction e + p — > e + p with large momentum transfer, the outgoing proton will 
be in a small-size configuration. But this small-size configuration is not detectable 
unless the outgoing proton has another particle to interact with. The role of the 
nucleus in this is that it serves as a laboratory to study the spatial configuration and 
spacetime evolution of a hadron produced in a hard exclusive reaction, by observing 
its interactions with the other nucleons. It is, however, important that the reaction 
be elastic or quasi-elastic; if the struck proton breaks apart, then it is not necessary 
for the outgoing particles to start in point-like configurations. The reason is that 
for a large-momentum transfer elastic reaction, all of the quarks involved have to 
be in a small region of space at the time of interaction in order for the quarks of a 



given particle to remain together after the collision (see Sec. 2.2.2 for more discussion 
of this). For an inelastic reaction, the outgoing quarks aren't required to remain 
together, and so they don't need to be close together at the time of the collision. 

Another feature of the large momentum transfer elastic reaction, e.g. e+p — > e+p, 
is that the produced small-size object (called an "ejectile") eventually must expand 
and become a normal-size proton, since a proton is what is detected. The small-size 
object is not an eigenstate of the strong Hamiltonian, and so evolves in time. For the 
reaction inside a nucleus, if the outgoing proton expands too quickly to its normal 
size, then it will experience normal-strength interactions with the other nucleons 
and no transparency will be observed. Therefore to observe Color Transparency, the 
outgoing proton should be moving fast enough that it has left the nucleus by the time 
significant expansion occurs. High velocity helps in two ways, the first of course being 
that it leaves the nucleus in a shorter time, and the second being that time-dilation 
slows down the rate of expansion as observed in the Lab compared to the rate in the 
proton's rest frame. 

Other reactions for which the same ideas hold include quasi-elastic proton scat- 
tering (p + A — > p + p + (A — 1)), pion photoproduction (7 + n — > n~ + p), pion 
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electroproduction, and the two which are explored in this thesis: pion elastic scatter- 
ing with proton knockout (tt+A — > n+p+(A— 1)), and vector meson electroproduction 
with proton knockout (7* + A — > V + p + (A — 1), V representing a vector meson, 
which in this thesis is the p). 

2.2.2 Formation of Pointlike Configuration 

The idea of why a large momentum-transfer elastic reaction involves a small size 
configuration of quarks can be seen heuristically from the example of the pion form 
factor, i.e. the reaction e + 7r — > e + tt, with large momentum transfer [2J. To lowest 
order in the electromagnetic interaction this occurs through emission by the electron 
of a virtual photon which is absorbed by one of the two valence quarks in the pion, 
the pion remaining intact. The 4-momentum of the virtual photon is denoted by q, 
with q 2 = —Q 2 . Loosely speaking, in the frame where the pion is moving fast, the two 
quarks are essentially moving parallel to each other at the same speed. When one of 
the quarks absorbs the large-momentum virtual photon (large Q) the direction of its 
motion is changed. If the pion is to remain intact, the momentum of the other quark 
must also change so that the two quarks are again moving in the same direction. 
This occurs by exchange of a gluon between the quarks. The 4-momentum-squared 
of this gluon is of order — Q 2 also, and hence the gluon is far off its mass-shell. By 
the uncertainty principle it can only exist for a time which is of order 1/Q. But in 
order to be absorbed by the other quark it must traverse the distance between the 
2 quarks, and so this distance must be less than order 1/Q. Thus in order for the 
pion to remain a pion, i.e. in order for the reaction e + tt — > e + tt to proceed, at the 
time the photon is absorbed the two quarks need to be closer together in space than 
a distance of order 1/Q. If they are then the reaction may proceed, and the outgoing 
"pion" will in fact consist of 2 valence quarks in a small-size configuration, with their 
separation being of order 1/Q. If the 2 quarks are farther apart than this, then they 
will separate from each other, with each quark eventually hadronizing, yielding an 
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inelastic reaction e + n — > e + X. Since the wavefunction of the quarks in the pion 
has some amplitude for the 2 quarks to be close together, the elastic reaction may 
occur. The outgoing ejectile then evolves over time, becoming the observed pion. 

Nonperturbative studies of realistic hadron models also show the formation of a 
small-sized configuration during a large- momentum transfer reaction [3]. 

2. 2. 3 Expansion from PL C 

Once a pointlike configuration is formed in a large-momentum transfer reaction, the 
system will expand until it reaches the "normal" size of the hadron; once it has 
reached its normal size the expansion ceases. In order to account for this, models 
must be used. The model used in the analysis in this thesis is called the "quantum 
diffusion model" [T^l |2]. In this model, the interaction cross-section of the outgoing 
object with the nucleons increases linearly with distance from the interaction point 
where the hard scatter occurred which produced the pointlike configuration (see Eq. 



3.3.27). This model is derived from perturbative QCD [5]: for a quark-antiquark 
system starting from a transverse size of zero, gluon exchange between the quark 
and antiquark proceeds until the system reaches the normal meson size. It is shown 
in [2] that the transverse area of the system (and hence its cross-section) increases 
linearly with distance traveled. The "naive" model of expansion would correspond 
to free quarks expanding from zero transverse size in both directions transverse to 
the momentum of the system. In this case the transverse area of the system would 
increase as the square of the distance traveled [T3]. Since the quantum diffusion model 
is derived from QCD (albeit perturbatively) , it is the model used in the calculations 
presented here. 

2.2.4 Experimental searches for CT 

The first dedicated experiment to search for effects of color transparency was in 1988 
at Brookhaven National Laboratory [7j. Quasi-elastic scattering of protons (i.e. the 



10 



reaction p + A —tp + p+(A — l), where a proton is knocked out of the nucleus by the 
incoming projectile) in various nuclei was observed, at incident proton momentum 
of from 6 to 12 GeV. The transparency, as a function of the 4-momentum transfer 
squared t, was observed to increase as \t\ increased, up to a point, but then the 
transparency decreased after that as \t\ was increased further. This behavior did 
not appear to agree with the predictions of color transparency, as the transparency 
should increase as \t\ is increased. However, their may be other factors at work in the 
elementary pp scattering cross-section, and several models were proposed to try to 
explain this behavior [1] . Another experiment was later performed [HI E] wherein the 
momenta of both outgoing protons was measured (in contrast to the first experiment 
where only one of the outgoing proton's momentum was analyzed). Similar results 
were obtained as in the earlier experiment, with the transparency first rising and then 
falling with \t\. 

In the (p, 2p) reactions, in order for a small-sized configuration to be formed it 
is necessary to have 6 quarks all localized in a small region, which may have a very 
small probability. The formation of a small-sized configuration may be more likely if 
fewer quarks are involved. Thus quasi-elastic electron scattering (e + A — > e + p + 
(A — 1)) may be a better candidate to observe color transparency. In this case, the 
elementary reaction e + p — > e+p is better understood also, being an electromagnetic 
interaction rather than a strong interaction. This experiment has been performed at 
SLAC [HI El H] with a range of momentum-transfer squared Q 2 from 1 to 8.1 GeV 2 . 
The results did not show any indication of color transparency. The observations 
agreed with the standard calculation which assumes that the outgoing object is a 
normal-sized proton with the usual free-space value of its cross-section of interaction 
with the other nucleons. 

There has been one experiment that can be said to show unambiguous evidence of 
color transparency. This was the diffractive dissociation into dijets of pions scattered 
from carbon and platinum nuclei [HI]. In this process a high-energy incident pion 
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strikes a nucleus, with the minimal quark configuration qq scattering coherently from 
the nucleus. The individual q and q then each form a jet of hadrons. Observation 
of two jets with large transverse momentum (transverse to the pion beam direction) 
indicates that the q and q had large relative transverse momentum and hence small 
transverse spatial separation. If the qq are in a pointlike configuration, then for 
forward scattering t ~ 0, since the scattering is coherent, the amplitude for scattering 
from a nucleus would be ~ A Ai, where M. is the amplitude for scattering from a 
single nucleon; hence the forward differential cross-section would depend on A as A 2 . 
In an experiment, what is measured is the integrated cross-section, a (A) = J ^-dt. 
For the coherent reaction, dadt is also proportional to the form factor of the nucleus, 
e tR\/3 w here r a j s the radius of the nucleus; this then gives an A-dependence of 
a(A) ~ A 4//3 [lj. This is to be contrasted with the expectation that a normal-size 
incident pion would undergo strong absorption from a nucleus, and so essentially only 
the nucleons on the surface would participate in the reaction; thus the A-dependence 
of the cross-section on a large nucleus would go like ~ A 2 / 3 . The result of the 
experiment pU] was a cross-section depending on A as A 1 - 55 , a clear indication of the 
effects of color transparency. 

One further experiment that has been performed, with somewhat inconclusive 
results, is the reaction 7 + n — > n~ + p in 4 He [TT| . The results show a momentum- 
transfer dependence that seems to indicate CT, with the transparency rising with \t\. 
However, better statistical precision is needed in order to be conclusive. 

Other candidate reactions are those involving production of vector mesons. These 
are discussed in the next chapter. 

As it would seem more likely to observe Color Transparency in reactions involving 
mesons, it would be of interest to measure the quasielastic scattering of pions from 
nuclei at large momentum transfer, i.e. tt + A — > ir + p + (A — 1). This reaction 
is the subject of this chapter. In the COMPASS experiment at CERN, pions with 
momenta of 200 GeV are produced. At this large momentum, the expansion of the 
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produced point-like configuration does not occur (due to time-dilation) before the 
pion escapes the nucleus. COMPASS should therefore be able to observe the effect of 
Color Transparency [12j. In order to calculate the cross-section, a formalism is needed 
which takes into account the initial- and final-state interactions of the pion, and the 
final-state interaction of the proton, with the spectator nucleons. As we are interested 
in high incident pion energy, the Glauber model provides such a formalism. In the 
following, the cross-section for the above reaction is calculated in the Glauber model. 
The result is the same as obtained in the usual Distored Wave Impulse Approximation, 
which has been used extensively to analyze proton knockout reactions [16] . The 
Glauber model can be easily extended to account for particle production processes, 
such as vector meson production, e.g. 7* + A — > V + p + (A — 1). In the following 
chapter this reaction is analyzed in the Glauber model. This represents a new result. 

2.3 Glauber model of high-energy hadron-nucleus scattering 

In order to calculate the scattering cross-section for a projectile incident on a nucleus, 
one needs a formalism which accounts for interaction of the projectile with more 
than one nucleon during its passage through the nucleus. The Glauber model of 
high-energy hadron-nucleus scattering [19] is a multiple-scattering model which is 
valid under certain conditions. The conditions are: 1) high energy of the incident 
particle, compared to the binding energy of the nucleons in the target nucleus; 2) small 
angle scattering of the projectile. Under these conditions the momentum transfer is 
mostly transverse, and so the longitudinal momentum transfer is neglected; the energy 
transfer from the projectile is also small, and so the energy transfer is neglected also. 
In order to calculate the scattering cross-section in the Glauber model, only knowledge 
of the free-space hadron-nucleon scattering amplitude and the wavefunctions of the 
target system is required. The Glauber model does not take into account the Fermi 
motion of the nucleons; for a projectile of high energy the Fermi motion should matter 
little. 
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The Glauber model takes advantage of the fact that high-energy elastic hadron- 
hadron scattering occurs at mostly small scattering angles; see e.g. [67] . For example, 
for proton-proton scattering at LAB incident momentum of 25 GeV, essentially all 
of the scattering events occur with \t\ < 1 GeV 2 , which corresponds to a center-of- 
mass scattering angle of 6°. This is also true for the usual eikonal approximation in 
potential scattering; the Glauber model is an extension of the eikonal approximation 
to include scattering from multiple scatterers in the target. In the Glauber picture, 
the nucleons' positions are fixed in place during the time that the projectile traverses 
the nucleus (the "frozen" approximation). Also, the projectile is assumed to scatter 
at most once from any individual nucleon. In between scattering events the projectile 
travels in a straight line. The Glauber result for the scattering amplitude is a sum of 
terms representing the possible multiple-scatterings of the projectile. The first term 
represents one elastic scatter, the second term represents two elastic scatters (from 
different nucleons), and so on up to a maximum of A elastic scatters, A being the 
nucleon number of the nucleus. 

The Glauber scattering amplitude is very similar to the Fraunhofer diffraction 
amplitude in optics, and can be interpreted in terms of diffraction [18]. For 2-body 
elastic scattering, with initial and final momenta k, k', we define the momentum 
transfer q = k — k', and for high-energy scattering we have q — q± where q_i_ is the 
component of q perpendicular to k. The scattering amplitude is then given by (see 
Appendix A) 

/(q) = — [ d 2 b e iq - b r(b) (2.3.1) 



2tt . 

which is the same as the expression for the scattering amplitude in Fraunhofer diffrac- 



tion, for scattering of an incident wave from an obstacle. In that case (see Fig. 2.1 ) the 
"profile function" T(b) is related to the amplitude transmitted through the obstacle 
(for incident amplitude 1) by t(b) = 1 — T(b). In the Glauber model, a wave incident 
on a system of n scatterers, each with profile function T(b), undergoes absorption and 



transmission through each scatterer (see Fig. 2.2). After it has passed through all of 
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i - r(b) 



(transmitted amplitude) 



b 



Figure 2.1: Diffraction of an incident wave by an obstacle. 



the scatterers, the transmitted wave has amplitude (at a given transverse position b) 



where Sj is the transverse position of the center of the j'th scatterer. The scatter- 
ing amplitude is then given by the 2-dimensional Fourier transform of 1 — £t t(b) = 
Ftot({b — Sj}). This is for scatterers at given positions. For a quantum system of 
scatterers (e.g. A nucleons in a nucleus) the scattering amplitude for the target sys- 
tem to remain in its initial state is given by the expectation value of the scattering 
amplitude in the initial state: 



n 



t tot (h) = n( i - r ( b - *)) 



(2.3.2) 



F 

- 1 it 




e iq - b r to4 (b, Sl ,...,s A ) (2.3.3) 




(2.3.4) 



For transition of the target system to a final state |/) the amplitude is 




(2.3.5) 



with the total profile function given by 



.4 



r tot (b,si 



-,BA)=i-n( i - r ( b - s i)) 



(2.3.6) 
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It is important to note that the states \i) and |/) are internal states of the target 



nucleus. (For a derivation of Eq. (2.3.5) for the case of a projectile scattering from a 



system of scatterers bound by a potential, starting from the Schrodinger equation for 
the projectile-target system, see [T§]). 




Amplitude of transmitted 
wave at this impact 
parameter b after passing 
through each obstacle: 



[i-r(b- Sl )] [i-r(b-. 



Figure 2.2: Diffraction from multiple scatterers. 



2.4 Quasi-elastic scattering: n + A — > n + p + (A — 1)* 

The amplitude for a transition from the target (A-nucleon) state |i) to state |/) is 
given in Glauber theory by 

Ffi = Y f d 2b ^' h [ d 3 r 1 ...d 3 r A <l>* f <l> i r tot (h,s 1 ,...,s A ) 5 3 (jJ2 r i) ( 2 ' 41 ) 
11 J J j=i 

where 

A 

r« = l-n[l-r(b- Si )] (2.4.2) 

and the profile function T(b) is related to the 7r-nucleon elastic scattering amplitude 
/(q) by /(q) = w~ J d 2 b e jq ' b T(b). Here q = k — k' with k = pion initial momentum 
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and k' = pion final momentum. And we take k along the positive z-axis, with b and 
Sj perpendicular to the z- axis, r,- = (sj, Zj) is the coordinate vector of the jth nucleon. 
$i(ri, r 2 , r^) and $/(ri,r 2 , r^) are the initial and final wavefunctions of the A- 
nucleon system; these are in principle exact eigenstates of the A-nucleon Hamiltonian. 
The final A-nucleon state |/) can be either a bound or continuum state. The delta 
function in the above equation enforces that the coordinates rj of the nucleons are 
relative coordinates; hence the target wavefunctions above are internal wavefunctions. 
For large A we can neglect the delta function. Note that for /(0) pure imaginary, the 
optical theorem in terms of T is J d 2 bT(h) = . 

For the case of proton knockout, the final state of the target nucleus consists of 
a continuum state, with one unbound proton. We will be using shell-model wave- 
functions for the initial target state. For the final A-nucleon state, we will assume a 
wavef unction of the form 

A (r!, r 2 , . . . , r A ) = x P (ri)<^_i( r 2, • • • , r A ) (2.4.3) 

where Xp( r i) is a scattering wavefunction for the proton of momentum p. Because of 
the (relatively) high energy of the outgoing proton, it is appropriate to use an eikonal 
wavefunction for the proton: 

X (r x ) = e ip ' ri e~^o 0dsa ^ tP( - ri+s ^ = e ip ' ri e~^ ap{ - ri \ (2.4.4) 

This represents scattering of the outgoing proton in the optical potential due to the 
other A — 1 nucleons. Therefore p in the exponential should be the nucleon density 
of the residual nucleus and hence would depend on the final state / of the residual 
nucleus. We will assume, however, that p is the same as the nucleon density of the 
initial nucleus, which should be valid for final states which are one-hole states or small 
excitations thereof. 

In an actual experiment, the outgoing pion and proton are detected, while the re- 
coiling residual nucleus is usually not detected. Therefore the cross-section of interest 



17 



is obtained by summing over all final states of the residual nucleus. In calculating the 
differential cross-section, however, the phase space factors depend in principle on / 
since the internal energy of the residual nucleus depends on /. For the case of a high- 
energy projectile, and proton knockout, it is legitimate to neglect dependence of the 



phase-space factors on /. Therefore we may simply square the amplitude Eq. (2.4.1), 
and then sum over all final states /, and we can use closure on the residual nucleus 

states: J2 f <f> A -i( r 2> • ■ ■ > r A)0A-i( r 2> • • • , t a) = llf=2 <5 (3) ( r j - r'j). The result gives a 
multiple scattering expansion, where the first term represents one hard scatter of the 
pion, with momentum transfer q, together with multiple soft re-scatterings as the 
pion travels through the nucleus; the second term represents two scatterings of the 
pion with momentum transfers qi and q 2 such that qx + q 2 = q, etc. For the case 
of quasi-elastic kinematics, where we have p ~ q, the outgoing proton has received 
almost all the momentum transferred from the pion. Therefore the higher-order terms 
(representing multiple hard scatters of the pion) should be negligible, and only the 
first term should be appreciable. Since in this term the pion only undergoes soft 
re-scatterings with the other A — 1 nucleons, and by assumption the outgoing proton 
also only undergoes soft re-scatterings (inherent in the eikonal form of the proton 
wavef unction), the final state of the residual nucleus should be a one- hole state of the 
target nucleus. And indeed this first term is identical to what is obtained if instead of 
summing over all final states of the residual nucleus we only sum over one-hole states. 
So let us evaluate that case. 

If we only sum over final states of the residual nucleus which are one-hole states, 
i.e. obtained from the initial A-nucleon state by deleting one single-particle state, 
then the initial wavefunction can be written (in the shell model) 

<f> A (ri, ...,r A ) = n (r 1 )0^_ 1 (r 2 , . . . , r A ) (2.4.5) 

with (fi n being a shell-model single-particle wavefunction, and the Glauber amplitude 
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is 



F (n) _ ik 

F ^ -to 



d 2 be^ h / d 3 n . . . rf 3 r AX ;(ri)0n(ri) ^fo, ...,r A ) T tot (h, {r,}) 



Separating out the terms in T tot which are independent of r 1; we have: 

A A 

r tot = 1 - J](i - r bj ) + r M J](i - r bj ) 

3=2 3=2 



(2.4.6) 



(2.4.7) 



= 1 - IJ(1 " r ^) + r e//(b, rx, ... , r A ) (2.4.8) 

3=2 

where 

r e// (b, n, . . . , r A ) = r 61 JI(i - r w ). (2.4.9) 

Because of the orthogonality of the single-particle wavefunctions \v an d 0«5 the terms 
in T tot that are independent of ri contribute zero to F/j. Hence only r e // contributes, 
and we have: 

F ^ = t J d2hel * h I ^•••^X;(ri)0 n (r 1 )|0Li(r 2 ,...,rA)| 2 r e// (b,{r J }) 

(2.4.10) 

For |0a-i| 2 , we assume an independent particle model, and write 

A 



A-ll r 2, • • • ,*A, 



3=2 



(2.4.11) 



with the single-particle density p\ normalized to 1. The nucleon density of the state 
4> A -\ is then p(r) = (A — l)pi(r). We then have 



/ 



d 3 r 2 . . . d 3 r A 



_ 1 (r 2 ,...,r A )| 2 r e// = r w f #r 2 ...tfr A Y\p 1 {T j ){l-T hj ) 

(2.4.12) 



3=2 



= r 



61 



J rf 3 r 2Pl (r 2 )(l-r(b-s 2 )) 



(2.4.13) 



= r 6 i g(b) 



(2.4.14) 
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and hence 



F { ; ] = £J d 2 be^ h J rf 3 r lX ;(ri)0 n (r 1 )r 61 g(b) (2.4.15) 

Since we are interested in large A, we can approximate g(h) by an exponential, 
as follows. We have 



J d 3 r 2Pl (r 2 )(l - T(b - s 2 )) = 1 - J rf 3 r 2Pl (r 2 )r(b - s 2 ) 

= 1 - J d 2 s 2 dz 2 p(r 2 )T(h - s 2 ) 



(2.4.16) 



Now the profile function T(b — s 2 ) is in general a sharply peaked function of its 
argument, peaked at s 2 = b, while in contrast the nucleon density p(r 2 ) is a much 
more slowly varying function. Hence we may approximate 

J d 2 s 2 dz 2 p(s 2 , z 2 )T (b-s 2 )~ J dz 2 p(h,z 2 ) J d 2 s 2 T{h - s 2 ) = ^/(0) T(b) 

(2.4.17) 

where T(b) = dzp(h, z) is called the "thickness function" [M], and /(0) is the 
forward pion-nucleon elastic scattering amplitude. For high-energy scattering, /(0) is 
almost pure imaginary, and so using the optical theorem we obtain 

J d 2 s 2 dz 2 p(s 2 ,z 2 )T(b - s 2 ) ~ X -a^ t T(b) (2.4.18) 

Hence 

r i ^-kN , X— 1 

g(b)~ l-^^T(b) ~ e -Kofnb) (2A1Q) 
J\ — 12 

Inserting this in the expression for Ffi, we have 

Fff = ^ J d 2 be^ h J d'nxl^Uri)^ T ^ (2.4.20) 

We may evaluate the integral over b at this point by using again the property of 
T(b — Si) that it is very sharply peaked at b = si, and write 

— / d 2 bT bl r(b) ^ e -\a^ T( Sl )ik I d 2 hT{h _ Si) {2A2l) 

2n J 2n J 

= e 4q ' Sl /(q) T (*i) (2.4.22) 
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and therefore 



Fff = /(q) / d'ne^xK^U^e-^ T(si) . (2.4.23) 

Finally, writing Xp( r i) = e~ ipri e~^ ap<yri \ we can write F{™' in terms of the missing 
momentum p m = p — q as 

= /(q) J rf 3 r 1 e- ipm - ri n (r 1 )e-5^( r i) e -| CT ^ T ( Sl ). (2.4.24) 



and 



\F { ; ] \ 2 = |/(q)| 2 y"d 3 r 1 d 3 r / 1 e-^-( ri - r i)0;(r' 1 )^ n (r 1 ) 



(2.4.25) 



Note that above we have used qsi ~ qri which is valid since q z <C |q_i_|- This result 
agrees with the usual Distorted Wave Impulse Approximation result for the amplitude 
of the (p, 2p) reaction [TB] if we identify the distortion factor for the incoming projectile 
as D(r) = e~2°"i*-°o dz p( s < z ) anc l the distortion factor for the outgoing projectile as 
D(r) = e~\ a $z dz'p(s,z')^ j s valid here since the scattering angle of the projectile 

is very small and so both integrals in the exponentials are along the same straight-line 
path. Note that these distortion factors are the same as one obtains in the eikonal 
approximation to the scattering wavefunction using an optical potential V op t(r) = 
— §erp(r) [65] . 

So now squaring and summing over all one- hole final states /, which is equivalent 
to summing over all occupied states n of the initial nucleus, we obtain 

A 

J2 \F$>\ 2 = |/(q)| 2 / ^ridV 1 e- ipm -< ri - r i> ^(r^i) 

n=l ^ 



(2.4.26) 

x e -|« P (ri)g-ia p (ri)g-|<7^T(si)g-|a^ T{s[) 



where 

1 A 

^O^S^WiN (2.4.27) 



A 

n=l 



is the shell-model one-body density matrix. 
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To evaluate this using shell-model wavefunctions, it's easiest to write it as 



El F /r } i 2 = i/(q)i 2 Ek (n) (p^ 



(2.4.28) 



n=l 



71=1 



where 
state a 



9 (n) (Pr. 



is called the distorted momentum distribution in the shell-model 



H6l. with 



9 (n) (Pr 



^ne" JPm ' ri 0n(ri)e~2 Q ^ r ^ e 



(2.4.29) 



For the proton knockout reaction, the transparency T is defined as the ratio of the 
measured 5-fold differential cross-section to the differential cross-section calculated in 
the Plane Wave Impulse Approximation (PWIA) jH El El [H] • This can be evaluated 
at a specific kinematic point, i.e. a particular value of the missing momentum p m , or 
it can be the ratio of the integrated cross- sections, integrated over some domain T> of 
p m . Thus 

da 

T(Pm) = (2.4.30) 



or 



T, 



Iv d3 P 



da 



m dE'dSl'dMr, 



V 



Sv d3 P 



da 



PWIA 



(2.4.31) 



m dE'dn'dn„ 



We call the latter the "integrated transparency" . At a given value of p m , the kinematic 



factors in the cross-sections cancel in the ratio Eq. |2.4.30| The 5-fold differential 
cross-section is proportional to Y^n=i l-^/^| 2 > anc ^ so we nave 

T{p m ) 



2^n=l \ r fi \PWIA 22n=l 



9^{Pm) 



2 

PWIA 



(2.4.32) 



The PWIA value of # (n) is obtained from Eq. ( |2.4.29[ ) by setting the attenuation 
factors equal to 1, which gives 

9pwiA{Pm) = J d 3 ne- iPm - ri Mn) = 0n(p m ) (2.4.33) 



22 



PWIA 



i.e. just the momentum space wavefunction of the nth state. Thus J2n=i 9^{Pm) 
is the momentum distribution of the initial nucleus. 

To incorporate effects of Color Transparency into our result, we note that the 



expression for the amplitude Eq. ( |2.4.24| can be interpreted as follows: the 

incoming pion strikes the proton in the nucleus at the position ri, which knocks the 
proton out; the proton suffers attenuation on its way out of the nucleus, beginning at 
the point r 1; as represented by the factor 



e -|ap(ri) 



(2.4.34) 



while the pion suffers attenuation on its way in (before the collision with the proton) 
from z = -oo up until the point z = z\ and on its way out (after the collision with 
the proton) starting at z = z\ until z = oo, as represented by 



e 2 



(2.4.35) 



Because the scattering angle of the outgoing pion is very small, it's legitimate to 
approximate its entire trajectory as being a straight line parallel to the z-axis. 

So now to include Color Transparency in the above result, we allow of£ and cxfj^ 
to depend on the distance of the given particle from the point where the pion struck 
the proton. Since the hard scatter occurs at the point r x = (si,2i), in the above 
formula we make the replacements 



.,„, , p{v 1 + s-p)ds -»■ / a^(s)p(r 1 + sp)ds 

'0 ^0 

/oo PZ POO 

dz'p( Sl ,z')^ dz'*l N (z')p( Sl ,z')+ dz'**J(z')p( Sl ,z') 
-oo J —oo Jz 

where the general form of the position dependent cr's is [T5] 



ttN 
'tot 



a 



tot 
WV 



0{h ~ z) 



z_ + n 2 (k 2 t ) ( ] 
h \t\ 



z 
h 



+ 9{z - l h ) 



(2.4.36) 
(2.4.37) 



(2.4.38) 
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where a to t is the free-space total cross-section (this model of the expansion from the 
pointlike configuration is called the "quantum diffusion model"). In this equation, n 
is the number of valence quarks of the hadron, while (k 2 ) 1 ^ 2 is the average transverse 
momentum of the quark in the hadron (taken to be (k 2 ) 1 ^ 2 = 0.35 GeV). Thus 
%^Cftjv is a measure of the transverse size of the hadron at the time of collision. The 



parameter lh (called the formation length; see Sec. 3.2.1 for more discussion of this) is 
the distance the hadron travels after the collision until it reaches its normal size. This 



i 



/, ~ p ■ c± M 2^ M 2 , where M n is the mass of a typical intermediate 



is estimated as L _ , . 
state n of the hadron [15] . In principle the quantities l n and l p can be different from 
each other, but since the relation lh — E \ Eh — M 2^ M 2 is only an estimate, we take 
here M 2 — = M 2 — M 2 = 0.7 GeV 2 for both l n and l p [TTj. 



Making the replacements Eqs. 2.4.36 and 2.4.37 in Eq. (2.4.29), we have 



#g(p m ) = yrf 3 r 1 e^ Pm - ri n (r 1 )e-5/o , ° CT r o TWp(ri+^) ds 



x e 2 



(2.4.39) 



For the case of tt + A — > tc + p + (A — 1)* with incident pion momentum |k| = 
200 GeV, the transparency T (Eq. (2.4.32)) was evaluated at p m = for a range 
of t from t = —1.5 to —10 GeV 2 , for the nuclei 12 C and 40 Ca. The result when the 



effects of Color Transparency are not included is obtained using Eq. (2.4.29) for g( 



in the numerator of T, while the result that includes effects of Color Transparency 



are obtained using Eq. (2.4.39) for gfW in the numerator. We shall call the former 



the "Glauber result" while the latter is the "CT result". The denominator of T is 
of course the same for both. The values of the free-space cross-sections used were 



a tot = 40 mb and = 25 mb which are valid for proton LAB momenta > 0.6 GeV 
and pion LAB momenta > 10 GeV |67j . 

For the wavefunctions <p n , harmonic oscillator wavefunctions were used. The os- 
cillator length b = */— was chosen so that the mean-square radius R 2 as calculated 



using the density from the wavefunctions, p(r) = \(f) n (r)\ 2 , was equal to the mean- 
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square radius B? as calculated using the Woods-Saxon form of the nuclear number 
density: 

p(r) = — ^ (2.4.40) 



r-R 

1 + e~ 



where R = 1.1 A 1 / 3 fm and a = 0.56 fm; po is determined by normalizing J d 3 rp(r) 
to A. The values obtained were b = 8.67 GeV- 1 for 12 C and b = 10.48 GeV' 1 for 
40 Ca. 
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Figure 2.3: Transparency T(p m ) for A = 12 and A = 40, for p m = 0. The bottom 
curves are the Glauber result; the top curves are the CT result. 



Fig. 2^ shows the calculated transparency at p m = 0. The effects of Color 
Transparency are very apparent. However, the value of the transparency is very 
sensitive to the value of p m . Fig. 2.4 shows the transparency, for A = 40, as a 
function of p mz for p mx = p my = 0, for p mz between and 400 MeV, for the Glauber 
case (i.e., not including effects of Color Transparency). 



2.5 Integrated transparency 



In the work of Benhar [T4], the integrated transparency Tx>, Eq. (2.4.31 ), is calculated 
by integrating over the entire region of p m such that the integrand is non-negligible. 
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Figure 2.4: Transparency for A = 40 as a function of p mz for p r 



P 



my 



0. 



As seen from Fig. 2.5, the distorted momentum distribution Y2 n =i ^"^(Pm)! 2 i s 



only significant for |p m | < 300 MeV. The PWIA value of this, which is the actual 
momentum distribution, is negligible for p m > 300 MeV, as the Fermi momentum 
of the nucleons in a nucleus does not much exceed this. The result in [H] for the 
integrated transparency, given without derivation, is 



Tv = \J d 2 sdzp(r)e- a ^e-^ T ^ 



(2.5.1) 



where the path of integration in ct p (r) is along direction of the outgoing proton's 
momentum p. In this section I derive this result, as well as the conditions under 
which it is valid. 



We are interested in calculating the integrated transparency, which is (Eq. (2.4.31 )): 



T, 



f v d 3 Pr , 



do 



dE'dn'dn^ 



g (n) (Pr, 



V 



Iv d3 Pm dE>Zm p J2n=l fv d ^m 9 {n) (p 



2 

PWIA 



(2.5.2) 



The second equality above is valid in the case where the phase-space factors in the 
differential cross-section are approximately constant over the domain D that is inte- 
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Figure 2.5: Numerator and denominator of T(p m ) versus p mz , for p mx = p my = 0. 
The dashed curve is the PWIA value (denominator) while the solid curve is the DWIA 
(numerator). These are both without inclusion of Color Transparency effects. 



grated over; in that case they factor out of the integral and cancel in the ratio. We 
wish to integrate over a region of p m corresponding to the entire range of momentum 
that the proton in the nucleus has. Therefore we integrate over all p m such that 
|Pm| < Vmax with p ma x — 300 MeV. For any |p m | larger than this, the integrand is 
negligible. Thus the numerator in the above equation is 



l num = j dZ v m j tfV^Me-^-^-^ <Pl{v\) n (n 



n=l 



(2.5.3) 

x e -§Mri) e -§Mr;) e -§<?nsi) e -KJ T(si) 



We have p = q + p m , and so if |q| ^> p ma x, then over the entire domain of integration 
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of p m we have p — q and we can re-write Eq. (2.5.3) as 

A 



71=1 ^ 

J d 3 p m e iPm < J rf 3 r 1 n (r 1 )e-5 ff ^ T ( Sl )e-^^ (ri) e- iPm ' ri 



x 

Pm ^Pmax 



(2.5.4) 

where we have set p = q in the a p 's. The last integral, over ri, is almost the 
distorted momentum distribution at p m for the state n; it differs from it in that it 
has ot q (ri) instead of a p (rx). We may therefore assume that this integral vanishes for 
Pm > Pmax] this is certainly true of the actual momentum distribution, which is just 
J d 3 ri<j) n (ri)e~ ipm ' ri . Therefore we may extend the upper limit on |p m | to infinity, 
with exact equality: 

n=l J (2.5.5) 
x J d 3 p m e iPm < J rf 3 r 1 n (r 1 )e-5^ T ( s i) e -i a ^ ri )e- iPm - ri 

where we now have integration over all p m G 3? 3 . Integrating over p m now gives a 
delta function, 5^'{tx — r^), and so finally we have, summing over n = 1, . . . , A 

I num ~ (2tt) 3 J d 3 r lP {r x )e-^ ^-O^M (2 5 6) 



Then the denominator of Eq. (3.4.1) is just the PWIA value of the above expresion, 



which is (27r) 3 J d 3 r\p(ri) = (2n) 3 A. Thus we have for the integrated transparency: 



Tv = jj rf 3 r 1 p(r 1 )e- ,J - T (»0 e -« 9 (n) (2.5.7) 



where the domain T> of missing momentum integrated over is all p m such that the 
distorted momentum distribution g( n \p m ) for all states n is non-zero (or at least 
non- negligible). 

It is important to note the essential assumption behind the preceding derivation: 
the momentum transfer |q| 3> p ma x where p max is the maximum momentum present 
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in the distorted momentum distribution. This allowed us to go from Eq. (2.5.3) to 



Eq. (2.5.4), which removed the dependence on p m from the a p 's. If |q| 3> p ma x does 
not hold, then it's certainly not the case that a p ~ a q . If, for example, |q| = p ma x, 
then in integrating over p m the direction of integration along the path of the outgoing 
proton in oe p (r) varies drastically. For r near the edge of the nucleus, that could make 
the difference between a p (r) being zero (if p = q + p m points radially outward) and 
Op(r) being significant (if p = q + p m points radially inward). Therefore in order for 



the expression Eq. (2.5.7) to be valid, we need to have the momentum transfer be 



much larger than the Fermi momentum of the nucleons in the nucleus. 



The result 2.5.7 is the result which is given in [15] as the semiclassical result for 
the transparency. We see that it does indeed have a semiclassical interpretation. 
The hard scatter with momentum transfer q occurs on a nucleon at the point (s, z), 
which knocks out the nucleon. The nucleon then propagates out of the nucleus. 
The incoming and outgoing pion, and the outgoing nucleon both suffer attenuation 
along their paths, which is given by the classical result for the attenuation of the 
intensity of a beam of particles passing through a material composed of scatterers of 
number density p(r). The position-dependent mean-free path of the particles in the 
material is then L(r) = (crp(r)) _1 where a is the cross-section of interaction, and so 
the attenuation factor starting from a given point r is just e~fo ds i L ( r + s p)] _ 



The integrated transparency Eq. (2.5.7) was calculated for A = 12, A = 40, and 



A = 208, for both the Glauber case (no CT effects included) and the CT case. In 



the CT case, the position dependent cross-sections Eq. |3.3.27| are used, as was done 
for the transparency T(p m = 0). The results are shown below in Fig. 2.6 It can 
be seen that the integrated transparency, for a given A and t, is smaller than the 
transparency T(p m = 0). In all cases the integrated transparency for the Glauber 
case is essentially independent of t, while for the CT case the transparency is much 
larger than for the Glauber case and increases markedly with \t\. 
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2.6 Concl usion 

We have calculated the transparency and integrated transparency for the proton 
knockout reaction it + A — > it + p + (A — 1)* within the Glauber theory, for the in- 
cident pion momentum of 200 GeV which is available at the COMPASS experiment. 
With the estimated values of the parameters that enter in the position-dependent 



cross-section Eq. 3.3.27, both the transparency T(p m = 0) and the integrated trans- 
parency T show large effects due to Color Transparency. In particular, for A = 208, 
even for modest values of \t\ the integrated transparency T is larger in the CT case 
than in the Glauber case by a factor of ~ 3 — 4, and increases substantially as \t\ 
increases. 
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Figure 2.6: Integrated transparency T. The bottom curves are the Glauber result; 
the top curves are the CT result. 
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Chapter 3 

COLOR TRANSPARENCY AND THE REACTION 

7* + A->p + P+(A-l)* 

3.1 Introduction 

In this chapter we calculate the transparency T(p m ) and integrated transparency T 
for the case of electroproduction of the p meson with proton knockout, 7* + A — > 
p + p + (A — 1)*. At the photon energies we are interested in, Glauber theory, 
modified to account for particle production, is valid, as it was for the case of pion 
scattering considered in the previous chapter. Electroproduction of the p provides 
another means of detecting the effects of Color Transparency In contrast to the 
purely elastic pion scattering considered in Ch. 2, for electroproduction there are more 
parameters that may be varied, namely the virtual photon energy v and virtuality 
Q 2 . These quantities, as well as a combination of them called the coherence length, 
l c = Q2 2 + m 2 ; can all affect the observed transparency. The coherence length plays 
an especially important role, since by varying its value the transparency T will vary 
even in the absence of any Color Transparency effects. Thus to observe an actual CT 
effect, one must keep the coherence length fixed. 



This chapter is organized as follows. In Sec. 3J2 electroproduction of vector 
mesons on a single nucleon is discussed. The Vector Meson Dominance model is 
introduced, and the coherence length and formation time for vector meson production 
are discussed. Existing experimental results in the search for CT involving vector 



mesons is discussed as well. In Sec. |3.3[ the Glauber formalism for particle production 
is presented. The amplitude for the reaction 7* +A — > p+p+ (A — 1)* is derived. Two 
limiting cases are analyzed, one for l c — > and one for l c — > 00, and it is shown that 
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for the case of l c — > oo the result reduces to the result for the pion elastic scattering 
case discussed in the previous chapter. The results for the values of the transparency 



T(p m = 0) are presented for several different values of A and l c . In Sec. 3.4, the 
integrated transparency T is calculated, for several different A and l c values. Sec. 3.5 
summarizes our results. 

3.2 Electroproduction of a vector meson on a single nucleon 

There are several pictures of electroproduction of vector mesons. In the Vector Meson 
Dominance model (VMD) [20J, the interaction of a real or virtual photon with a 
nucleon proceeds with the photon first fluctuating into a (virtual) neutral vector 
meson (i.e. a meson with the same quantum numbers as the photon), followed by the 
virtual vector meson scattering elastically from the nucleon. The elastic scattering of 
the virtual meson on the nucleon puts the meson on its mass shell. The amplitude 
for the production process 7* + iV — > N + V is then proportional to the elastic 
scattering amplitude for V + iV — > V + N. In this picture, the physical photon is a 
superposition of a bare photon state and vector meson states (the bare photon state 
would be the real photon state in the absence of the strong interaction). Thus at 
a given photon energy and q 2 , and a given momentum transfer t = (p 7 — p p ) 2 , the 
production amplitude is simply proportional to the elastic scattering amplitude. 

Electroproduction of vector mesons can also be described in terms of quarks, using 
QCD. The virtual photon fluctuates into a virtual qq pair, which propagates over 
a distance called the coherence length (determined by the energy-time uncertainty 
principle) before scattering elastically from the nucleon, which puts the qq pair on 
the mass-shell of the vector meson. The qq state then evolves over time to form the 
final real vector meson state. The transverse size of the qq that is produced by the 
virtual photon goes as r_|_ ~ 1/Q [1], so the larger Q is, the smaller is the size of the 
produced qq. In the limit of Q — > 00 the size goes to zero: a point-like configuration. 
Thus for large Q 2 the produced object should have vanishing interactions with the 



33 




Figure 3.1: Vector meson production in the VMD model. The incoming photon 
fluctuates into a virtual vector meson (V) which then scatters elastically from the 
nucleon. 



other nucleons and the transparency should approach 1. 
3.2.1 Coherence length and formation time 

There are two length scales (or time scales) of relevance to vector meson production, 



the coherence length and the formation time (see Fig. 3.2). The distance that the 
virtual hadronic fluctuation of the photon can travel in the LAB frame (target nucleon 
or nucleus at rest) is known as the coherence length pp. The energy-time uncertainty 
relation is used to determine this distance. For a photon and virtual meson with 
the same momentum k, the difference in energy between the photon and the virtual 
meson is 

A E = yjk* + ml - v = yjk* + ml - ~ Q " ~ Q " (3.2.1) 

where we've assumed v ^> Q, my. For this high-energy case, the velocity of the 
vector meson is essentially c, and so the energy-time uncertainty relation gives the 
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coherence length as 




(3.2.2) 



nucleon 

Figure 3.2: Coherence length (l c ) and formation length (If) for vector meson produc- 
tion. The incoming photon dissociates into a qq pair which then interact with the 
nucleon by exchanging gluons. 

For vector meson production in a nucleus, while the virtual hadron or qq is prop- 
agating over the distance l c it may interact with nucleons and be absorbed, before it 
has a chance to undergo the elastic scatter which puts it on mass-shell. These Initial 
State Interactions (ISI) therefore affect the measured production cross-section in the 
nucleus. In general, as l c increases, the probability of absorption increases and so 
the measured production cross-section in a given nucleus should decrease. Thus the 
production cross-section at low energy (small v) should be larger than the production 
cross-section at high energy (large u), for a given Q 2 . Or conversely, for a given z/, as 
Q 2 is increased, l c will decrease and therefore the measured production cross-section 
should increase. This effect mimics the effect of Color Transparency. Therefore in 
order to detect effects of CT, the coherence length should be kept fixed in a given 
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experiment. 

The formation time is the time scale over which the virtual meson or qq pair 
develops into the final real vector meson state, after scattering from the nucleon. The 
scattering with the nucleon puts the virtual meson or qq pair onto the mass shell of 
the vector meson. At the time of scattering the transverse size of the qq is small, and 
as it propagates away it evolves into the final meson state. This time can be estimated 
by considering the on-mass-shell small-size qq pair as a superposition of hadron states, 
namely the final real vector meson state and the next higher-mass meson state pQ. 
Then the energy-time uncertainty principle in the rest frame of the outgoing meson 
gives 

At = , (3.2.3) 

myi — my 

while in the LAB frame this is time-dilated by the factor 7 = —, — — r, where 

J 1 ^(rn v /+rn v ) 7 

^(my + my) is the average mass of the two states, and so the formation time or 
length in the LAB (assuming (3 ~ c) is 

lh = 2 2PV 2 - (3.2.4) 
m v , - m v 

3.2.2 Experimental results for p electroproduction in nuclei 

There have been several searches for evidence of Color Transparency in electropro- 
duction of p mesons in nuclei. At Fermilab in 1995 [2T], high energy muons were 
scattered from nuclei to produce p's. It was thought that CT was observed because 
the transparency, for a given A, increased as Q 2 was increased. However, in this 
experiment the coherence length was not held constant as Q 2 was increased, so it 
is difficult to draw conclusions from their data. A later experiment at DESY was 
conducted to explicitly measure the coherence length effect [22]. It was observed, as 
expected, that the transparency decreased as l c was increased, in p electroproduction 
in 14 N. The Q 2 values for this experiment were such that no CT effects should occur, 
i.e. the produced object would interact with the full p-nucleon cross-section. Hence 
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any dependence of the transparency on l c was not an indication of CT. This was a 
clear indication that any attempt to detect CT in vector meson electroproduction 
must look for effects while holding l c constant. Another experiment at DESY [23] 
was performed, where the transparency as a function of Q 2 was measured for differ- 
ent values of l c . There appeared to be an increase in the transparency as Q 2 increased, 
although the number of events at each fixed value of l c was not large, and so bet- 
ter statistics are needed. Finally, the most recent experiment to search for CT in p 
production was at JLAB [21]. In this experiment, the coherence length varied from 
~ 0.5 fm to ~ 0.85 fm. For this range of coherence length, the qq is produced 
essentially right at the location of the nucleon that it scatters from, and so there are 
no Initial State Interactions. The transparencies on 12 C and 56 Fe were measured for 
Q 2 from 1.0 to 2.3 GeV 2 . The transparencies appeared to show an increase with Q 2 , 
although statistics again were low. 

3.3 The Glauber formalism for particle production 

For the case of particle production, the profile operator T to t now depends on lon- 
gitudinal momentum transfer [25J. The reason is that for forward production on a 
single nucleon there is necessarily non-zero longitudinal momentum transfer due to 
the difference in mass between the incident particle and the outgoing particle. For 
the case of 7* + N — > p + iV at high energy, the energy transfer from the photon to 
the nucleon can be neglected, and so for an incident photon of momentum k, energy 
v 1 4-momentum squared k 2 = v 2 — k 2 = — Q 2 , and an outgoing particle of mass M, 
energy u, and momentum k' with momentum parallel to k (i.e. forward production) 
conservation of energy gives 

\k'\ 2 + M 2 = |k| 2 -Q 2 . (3.3.1) 
With the longitudinal momentum transfer = |k| — |k'|, we have 

M 2 + Q 2 = |k| 2 - |kf = 2|k| q L + 0(q 2 L ) (3.3.2) 
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and so for v ^> Q we have 

Q 2 + M 2 

q L = . 3.3.3 

Iv 

This longitudinal momentum transfer modifies the profile function T(b) [64], due 
to the phase difference between the incident (photon) wave and the outgoing (meson) 
wave. Consider vector meson production on a nucleon located at z 3 - (with the incident 
photon along the z-direction) . The phase of the transmitted wave at a point z > Zj 
equals the phase of the incident (photon) wave at z 3 - plus the change in phase of the 
transmitted wave as it propagates from Zj to z. Thus the transmitted wave at the point 
z is e lkz j e lk '{ z - z j) — e ik'z e tzj(k~-k') _ e ikz e iq h zj p or e i as ^j c scattering of a projectile, the 

wave at point z would just be e lk ' z . Therefore the phase difference of the incident and 
transmitted waves is just qLZj, and so the profile function for production on a nucleon 
at (sj,Zj) is T 7 (b — 8j)e tqLZj = Tj J e iqLZj . (The notation here is the same as in the 
previous chapter: b and Sj are two-dimensional vectors perpendicular to the incident 
photon's momentum direction, and z and Zj are coordinates along the z-axis which 
is parallel to the incident photon's momentum). Here T 7 is related to the production 
amplitude f lV (q) for 7* + iV — > V + N (where q is the transverse momentum transfer) 
by 

/ 7 ^(q) = — [ d 2 be l(l h T 7 (b) (3.3.4) 
2n J 

which is the same as for the elastic scattering case. So we have also 

T 7 (b) = J d 2 qe-^ h r v (q), (3.3.5) 

giving T 7 in terms of f jV . 

Taking into account qL, the total profile operator T tot represents production of the 
vector meson on a nucleon at (sj, zj), followed by any number of re-scatterings of the 
produced meson on the other nucleons. Hence T tot has the form [OH EE] : 

A 

T tot (b, r 1 ,...,r A ) = J2 r l/ qLZj Hi 1 - T bm mj \ (3.3.6) 

j=l m^j 
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where T bm = T(b — s m ) is the profile function for elastic meson-nucleon scattering, 
and 9 m j = 6{z m — Zj) ensures that any elastic scattering of the produced vector meson 
occurs after the meson has been produced (for high-energy scattering, the waves are all 
"moving forward" , which is along the z-direction, and so "later in time" is equivalent 
to "farther along in the z-direction") . 

As in the pion case (Ch. 2), we will sum over the residual nucleus final states 
which are one-hole states of the initial nucleus, and so the scattering amplitude is 
again 

Ff = ^ J d 2 be^ h J c?r 1 ...^;(r 1 )0 n (r 1 )|^_ 1 (r 2 ,...,r A )|V«(b,{rJ) 

(3.3.7) 

Expanding out T tot into terms that depend on ri and terms that don't, we have 



v tot = J] [i - r bk e kl ] - r bl £ r;,/"' W'i, JJ [1 - v hk e kj \ (3.3.? 

3=2 k^l, j 



k=2 



+ Y, T V 9LZ3 II [i-r»^] 

3=2 Ml, j 



(3.3.9) 



and note that the third term is independent of r x and so contributes zero to Ffi due 
to orthogonality of Xp an d <f> n - Therefore we have 



F { ; ] = ^ / d%e^ I ^...^^(^^(rO^Li^,...,^)!^^^}) 



(3.3.10) 



where 



A A 

= ^ i9lzi ui 1 - t *m - r ^ e r ','"' 1 n t 1 - ( 3 - 3 - n ) 

3=2 Ml, 3 



k=2 



Taking again an independent particle model for the residual nucleus, so that 



T A ) =Y[f>l{Tj), 

3=2 



(3.3.12) 
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the second term in Eq. (3.3.11 ) contributes A — 1 equal terms to Fu. Performing the 



integral over y 2 , . . . , r A we obtain 



d 3 r 2 . . . d 3 r A 



y A _ 1 {r 2 ,...,r A ) 



r e// (b,{ rj -}) 



Y^e^g^h) - r 61 / d 3 r 2 p(v 2 )Tj 2 e^ 6 12 g 2 (b) 



where 



and 



gi(h) = 1— / d 2 s I dzpi(s, z)T(b — s 



g 2 (b) = 1 - / d 2 s / ^pi(s,z)r(b - s 



Z2 



A-2 



In the large- A limit, we have 



g 2 (b)~ e -^ T ^ 



(3.3.13) 
(3.3.14) 

(3.3.15) 

(3.3.16) 

(3.3.17) 
(3.3.18) 



where Tj(b) = J 2 °° dz'p(b, z') is the "partial thickness function". 

As in the pion case, we can utilize the fact that the profile functions are sharply 
peaked and the other factors are relatively slowly varying. For a slowly varying 
function /(r) we thus have, to good approximation, 

J d 2 sf(s, z)T(s - a) ~ /(a, z) J d 2 sT{s - a) (3.3.19) 

and similarly for r 7 (s — a). 

Using the above approximation, we may integrate over s 2 in Eq. 3.3.14, using 

/2ir 
d 2 s 2 V(b ~s 2 ) = - 



(3.3.20) 



and then we may integrate over b in Eq. 3.3.10 by setting b = Si everywhere except 
in the profile functions and e iqb , with the result: 



J d 2 s 1 dz 1 e- ipm - ri e-^ ap{ri) (l) ri (v 1 ) 

2fc J-oo 



i} /(q; 



(3.3.21) 
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(a) First term of Eq. |3.3.21 



(b) Second term of Eq. |3.3.21 



Figure 3.3: Pictorial representation of the two terms in the amplitude of Eq. |3.3.21 



Here we have written the result for in terms of the missing momentum p TO , which 
is defined by 

p m = p - k + k' = p ± - q + (p z - q L )z (3.3.22) 
where p is the momentum of the outgoing proton. 



The physical interpretation of the two terms in 3.3.21 is as follows. The first term 
in parentheses corresponds to production of the vector meson on nucleon 1 at position 
(si, z±) with transverse momentum transfer q, nucleon 1 being therefore knocked out. 
The second term in parentheses corresponds to forward production of the vector 
meson on nucleon 2 at position (s^,^); the produced meson then propagates in the 
^-direction until the point (si,zi) where it scatters elastically from nucleon 1 with 
transverse momentum transfer q to nucleon 1, nucleon 1 being knocked out. In both 
cases the vector meson suffers attenuation beginning at the point where it is created 
as a physical meson through interaction with a nucleon (either at (si , Z\) for the first 
term or at (si, z%) for the second term), while the proton suffers attenuation beginning 
at the point = (sj, Z\) where it was located when the vector meson struck it. The 
total amplitude is the sum of these two amplitudes; hence the square of the amplitude 
contains interference between the two amplitudes. 
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The result Eq. 3.3.21 is the Glauber theory result for the scattering amplitude for 
7* + A — > p + p+ (A — 1)*, for the case where the final state of the residual nucleus is a 
one-hole state of the initial nucleus. To obtain the differential cross-section, summed 
over all one-hole states, we would square i* 1 !™ multiply by the appropriate phase- 
space and flux factors, and then sum over n = 1 to A. For the high-energy case we 
are considering, we may consider the energies of the outgoing particles to be essentially 
independent of n. In that case, the phase-space and flux factors are independent of 
n, and so we may just sum \Fj? | 2 over n. No inclusion of Color Transparency effects 
has been made up to this point, since the cross-sections a tot that appear in it are the 
measured free-space cross-sections, and the elastic p-nucleon rescattering amplitude 
/(q) that appears is also the free-space elastic amplitude. Hence the outgoing meson 
or proton interacts with the other nucleons with the full free-space interaction cross- 
section. To include CT effects, these cross-sections, and the elastic amplitude /(q), 
must be modified to account for the smaller size of the outgoing hadrons compared 
to their usual sizes. 

The result for Y2 n l-^/i | 2 > where n is summed only over one-hole states, is identical 
to the result one would obtain if instead one summed over all final states of the residual 
nucleus (the incoherent cross-section) but only kept the terms corresponding to a 
single rescattering of the produced vector meson on a proton, and neglected terms 
where the vector meson rescatters two or more times on different nucleons. The 
experimental situation, wherein the recoiling nucleus is not detected, corresponds to 
summing over all final states of the residual nucleus. However, because of the exclusive 
nature of the reaction, if p m is small, then the outgoing proton's momentum p — q and 
so only a single rescattering of the p can have occurred, where the entire momentum 
transfer q was delivered to the detected proton. Multiple rescattering terms in this 
case should be negligible, and so we need only sum l-Fj™^ 2 over one-hole final states. 



This implies that the transparency T using the result Eq. 3.3.21, which neglects any 
Color Transparency effects, will show very little dependence on the 4-momentum- 
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transfer-squared t ~ — q 2 . Any significant variation of T with t will be due to Color 
Transparency. 

Two limiting cases of the above result are of interest. For the case of large qi, 
which corresponds to a small coherence length l c , the above result simplifies. Taking 



qi — > oo, the second term in Eq. (3.3.21 ) is zero because of the oscillating exponential. 
So in that case, 



fi 



r v (q) I cfMzie- iPm - ri e- 



5«p( r l), 



>n(ri)e 



Ti(si) 



(3.3.23) 



This is similar to the result for the pion case, Eq. (2.4.24), except that the vector 



meson only undergoes attenuation (the factor e~?°*°* Tl ( Sl )) starting at the point z\, 
which is also the point where the initial proton was when it got knocked out. There 
is no attenuation before this point; this agrees with the small coherence length, which 
means that the photon fluctuates into the vector meson essentially at the same point 
where it interacts with the proton with momentum transfer q. For the case of the pion 
(Ch. 2), the incoming pion can of course interact all along its incoming trajectory, 
and hence its attenuation factor e~\.^ includes integration from z = — oo to 
z = oo. 

The other limiting case of interest is for q^ = 0. In this case the profile operator 



T t ot (Eq. (3.3.6)) is equal to the T tot for the pion case, Eq. (2.3.6), if we take T 7 = T. 
This is easily shown for A = 2 (and then proved for arbitrary A by induction on A). 
For A = 2 we have: 



r tot (b, n, r 2 ) = r 61 [i - r 62 # 21 ] + r b2 [1 - r bl e 12 ] 
= r^x + v b2 — r fel r fe2 #2i — r b2 r ;)1 ^ 1 2 
= r M + r b2 - Y bl v b2 = i - (i - r 61 )(i - v b2 ) y/ 



(3.3.24) 



Therefore the result Eq. (3.3.21) should also reduce to the result for the pion case, 
Eq. (2.4.24), when we set qi = and T 7 = T, and indeed it does: for T 7 = T we 



have / 7V (q) = /(q), and for qi = the second term in parentheses in Eq. (3.3.21) 
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becomes 



1 VNr/ \ f Zl r t \ -\aYJ f°° dz' p(si,z') 

-7ftotf{<\) / dz 2P {s 1 ,z 2 )e 2 tot ^ 



x ^°~YJt ' dz 2 



I V/,,,., / .,. / - ^ f 2 .... 



to* 

= /(q )(_ e -Kf T l(Sl ) + e -|<f T( Sl)) _ 

Note that the optical theorem was used to relate f lV (0) = /(0) to cr^. Thus we 
have 



(3.3.25) 



F (n) = flV ^ I d 2 Sl d Zl e- iPm - ri e-^ ri) ' 



in agreement with the pion result. 



(3.3.26) 



3. 3. 1 Inclusion of Color Transparency effects 



Effects of Color Transparency can be incorporated into the result Eq. (3.3.21) by 



including position dependent cross-sections from the quantum diffusion model [15112]. 
In this model, the total cross-section of interaction of the outgoing hadrons with a 
nucleon in the nucleus is 



a 



eff, 
hN 1 



Z,t) 



a 



tot 
hN 



0(k - z) 



+ 



h 



(3.3.27) 



Here z is the distance the hadron has traveled from the point where the hard hadron- 



nucleon interaction (with 4-momentum-transfer-squared t) occurred (Fig. 3.4), a 



.tot 
'hN 



is the free-space total hadron-nucleon cross-section, n is the number of valence quarks 
of the hadron, and (k 2 ) 1 ^ 2 is the average transverse momentum of the quark in the 
hadron (taken to be (k 2 ) 1 ^ 2 = 0.35 GeV). Thus ^p-cr^y is a measure of the transverse 
size of the hadron at the time of collision. The parameter 1^ (the formation length) is 
the distance the hadron travels after the collision until it reaches its normal size. This 
is estimated as h, — F 1 F — , Z v \ t i , where M n is the mass of a typical intermediate 



En — Eh 



Ml-Ml ' 
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state n of the hadron [IS]. In principle the quantity lh can be different for the pion 
and the proton, but since the relation lh — E \ Eh — f i is on ly an estimate, we 



M 2 -M? 



take here M 2 — = M 2 — Ml = 0.7 GeV 2 for both l n and l p [H]. 




Nucleus 



Hard scatter 

Figure 3.4: Formation length (lh) for vector meson production, z is the distance of 
the outgoing hadron from the point where the hard scattering occurred. 



The expression Eq. 3.3.27 is used for the cross-sections that appear in the expo- 



nentials in Eq. 3.3.21 The amplitudes / 7l/ (q) and / 7,/ (0) that appear in Eq. 3.3.21 
are the same as the measured free-space production amplitudes. However, the elastic 
rescattering amplitude /(q) must be modified to include the effects of Color Trans- 
parency. For large enough Q 2 , the qq pair produced at the point (si,Z2) will be in 
a pointlike configuration; it will then expand as it propagates, and scatter elastically 
from a nucleon at z 2 ; if z 2 is close enough to Z\, the scattering amplitude /(q) of the 
qq pair on the nucleon will be smaller than that of a normal p meson. Therefore the 



scattering amplitude /(q) in Eq. 3.3.21 should be replaced by [28 



f PLC (z 1 -z 2 , q ,Q 2 ) = f( q ) 



a e v f i( Zl -z 2 ,Q 2 ) G v{t 



T tot 
'VN 



a 



tot 
VN 



G v (t) 



(3.3.28) 
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where Gy(t) is the p-meson form factor, and t ~ — q 2 , and /(q) is the measured free- 
space elastic p-nucleon scattering amplitude. This form for f PLC is derived using the 
optical theorem (and assuming /(0) is pure imaginary) together with the empirical 
result EE] that the differential cross-section for hadron-nucleon scattering satisfies 

j hN-^hN 

^ Gl(t)G%(t). (3.3.29) 

in terms of the form factors of the h and N. 

Thus the result for the scattering amplitude including Color Transparency effects 

is 

= J rf 2 s 1 ^ 1 e- iPm - ri e-^p( r i)0 n (r 1 ) 
x (/^(q)e-W(^i) - ^P V (0) p dz 2 p( Sl ,z 2 ) e ^~^ e~^ v ^ f PLC {z u z 2 , q, Q 2 



(3.3.30) 



where 



roc 

ttp(ri) = / a e p f N f (s,t)p(r 1 + sp)ds (3.3.31) 
Jo 

poo 

a v {s l ,z 1 )= dz'a e v f l{z' - z 1 ,t)p{s 1 ,z') (3.3.32) 

J 21 
POO 

a v (s 1 ,z 2 )= dz'a e Jl(z' - z 2l Q 2 )p(s 1} z'). (3.3.33) 

J Z2 

These expressions for ay reflect the fact that the transverse size of the initial qq (at 
z 2 ) is determined by 1/Q 2 , while the transverse size of the outgoing qq and proton, 
after the hard scatter from the proton at (si, z±), is determined by l/\t\. 

The transparency T was calculated for 12 C and 40 C a at p m = for kinematics 
corresponding roughly to those in the JLAB proposal for electroproduction of p in 
nuclei [SB]- The same harmonic oscillator nuclear wavefunctions were used as were 
used for the pion scattering case in Ch. 2. The free-space cross-sections used were 
o*% = 40 mb, and cr£f = 25 mb P]. 



Graphs of T(p m = 0) vs. Q 2 are shown in Figs. 3.5 and 3.6 It is important 
to note that the transparency as a function of t is calculated for fixed v and Q 2 , so 
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Figure 3.5: Transparency T(p m ) for p m = for A = 12, t — — 2 GeV 2 , l c — 5 fm . 
The bottom curves (gray) are the Glauber result; the top curves (black) are the CT 
result. The value of the elastic p-nucleon t-slope parameter b used in the calculation 
is indicated for each graph; VMD corresponds to 6 7 y = b. 
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Figure 3.6: Transparency T(p m ) for p m = for A = 40, t = -2 GeV 2 , l c = 5 fm . 
The bottom curves (gray) are the Glauber result; the top curves (black) are the CT 
result. The value of the elastic p-nucleon t-slope parameter b used in the calculation 
is indicated for each graph; VMD corresponds to 6 7 y = b. 
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that the coherence length l c is held constant. If the coherence length varied, this 
could mimic Color Transparency because as l c gets smaller the attenuation due the 
Initial State Interaction of the vector meson (before the hard scatter) decreases since 
the vector meson propagates a smaller distance before undergoing the hard scatter; 
this would cause the value of the transparency T to increase as l c decreases. The 



production and elastic scattering amplitudes / 7y (q) and /(q) in Eq. (3.3.30) and 
Eq. J3.3.28| ) were taken to be of the form p v (q) = A lV e^ vt and p v (q) = Ae^ ht 
(where t = — q 2 ) with the parameters A^y, b 7 y, A and b taken from experimental 
data. The t-slope b for elastic p-nucleon scattering has been measured to be between 
7 and 8 GeV~ 2 [69J, while the t-slope for the production amplitude varies with Q 2 . 
The available electroproduction data [2Z] are at higher virtual photon energies than 
are considered in this paper, but the values of & 7 y(Q 2 ) measured in that experiment 
were what were used in our calculations. Calculations were done for b = 7 GeV -2 
and for 6 = 8 GeV -2 with 6 7 y depending on Q 2 . For comparison, calculations were 
also done assuming the validity of Vector Meson Dominance, in which case 5 7 y = b 



and the transparency T(p m ), Eq. (2.4.32), is independent of the value of b since 
both numerator and denominator are proportional to e bt . The expected properties 



of the transparency are evident in Figs. 3.5 and 3.6 For a given value of Q , the 
transparency (both Glauber and CT results) decreases with increasing A. For a given 
A, as Q 2 increases the transparency in the CT case increases, which is also expected. 
However, for the Glauber case, the behavior of T as Q 2 varies is more sensitive to the 
values of b and 6 7 y that are used. Some of the dependence of T on Q 2 is also due to 



the dependence of ot p (r) on kinematics through the relation Eq. (3.3.22). 
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3.4 Integrated transparency 

As in Sec. 2.4.2, the experimental situation corresponds to a range of values of the 
missing momentum p m . The integrated transparency is again 

2 

(3.4.1) 



m dE'dn'dn p 

dapwiA 



V 



En=l L d3 P- 



p. 



fyd 3 P . 

Iv d3 Pm 2§mt p Iv d'Pm FW ( Pm ) 



2 

PWIA 



Following the same steps as in the case of pion scattering, Sec. 2.4.2, if we integrate 
over p m up to p ma x — 300 MeV, we may set p = q in a p ; then assuming that the 
momentum distribution is zero for p m > p m ax, we may extend the upper limit to 
infinity, p m ax oo. For the denominator we obtain simply (2tt) 3 A \f lV (q)\ 2 . For 
the numerator we obtain 3 terms: 

(27r) 3 |/ 7V '(q)| 2 / dh^pir^e-^ih^) + ^(n) + /i 3 (n)) (3.4.2) 

where 



h^n) = e - av{si ' zi) 



(3.4.3) 



4tt / 7 (q)/ 7 {0) „-i av ( sl ,z 1} I ' dz 2 p(s l ,z 2 )e- 1 2 av ^' z ^cosq L (z 1 -z 2 )f PLC (zi,z 2 ,q) 



h(ri) = (- 



2ir\* \p v (0)\ 



% f Zl dz 2 r dz 3 p(s 1 ,z 2 )p(s 1 ,z 3 )e-^ s ^e- ± ^ s ^ C osq L (z 2 -z 3 ) 

q)l J-™ J-™ 



(3.4.4) 

-\a v {s 1 ,z 2 ) „-=ay(si,2 3 ) 



xf PLC (z 1 ,z 2 ,q)r^(z 1 ,z 3 ,q 



*PLCi 



Thus we have for the integrated transparency 

Tv = \\ rf2 ^iP( r i) e ~ Mri) (Mri) + h 2 {r x ) + h z {r x )). 



(3.4.5) 



(3.4.6) 



This simplifies somewhat if we assume the validity of Vector Meson Dominance for the 
relation between the free-space production amplitude p v (q) and the free-space elastic 



scattering amplitude /(q) (which appears inside f PLC ; see Eq.3.3.28). Assuming that 
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the high-energy amplitudes are purely imaginary, use of the optical theorem then 
gives: 



h 2 (rx) 



a VN e -§ay(»i,*i) 



G v (t) 
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A\G v (t) 
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(3.4.7) 
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where 



h(z) = 



a 



tot 
VN 



0{k - z) 



g n2 ( k t) (, _ i 
h Q 2 v / 



+ 0fc - 



(3.4.9) 



The form factor Gy used in evaluating Eq. |3.4.6| was taken to be the same form factor 
as for the pion: 



G v (t) 



1 



1 - J/0.59' 



(3.4.10) 



for t in GeV 2 . 



The 3 terms of Eq. 3.4.2 or Eq. 3.4.6 are represented pictorially by the same 



diagrams as in Fig. 3.3 The term with hi is the square of the diagram in Fig. 3.3(a 



and represents incoherent production from nucleon 1; the term with h 2 represents 



interference between the diagrams of Fig. 3.3(a) and Fig. 3.3(b) with interference 
between production on nucleon 1 and nucleon 2; and the term with h 3 is the square 



of the diagram in Fig. 3.3(b) , which represents interference between production on 



nucleon 2 and production on a different nucleon 3, with incoherent scattering from 
nucleon 1. 

The integrated transparency was calculated for A = 12 and A = 40, for a range 



of values of t and Q . In Figs. 3.7 - 3.10 the transparency is shown for fixed t as a 
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(c) b = 8 GeV~ 2 



Figure 3.7: Integrated transparency T for A — 12, t — —2 GeV 2 , l c — 2 fm. The 
bottom curves (gray) are the Glauber result; the top curves (black) are the CT re- 
sult. The value of the elastic p-nucleon t-slope parameter b used in the calculation is 
indicated for each graph; VMD corresponds to b 7 v = b. 
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(c) b = 8 GeV" 2 



Figure 3.8: Integrated transparency T for A — 12, t — — 2 GeV 2 , l c — 5 fm. The 
bottom curves (gray) are the Glauber result; the top curves (black) are the CT re- 
sult. The value of the elastic p-nucleon t-slope parameter b used in the calculation is 
indicated for each graph; VMD corresponds to b y v = b. 
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Figure 3.9: Integrated transparency T for A = 40, t = —2 GeV 2 , l c — 2 fm. The 
bottom curves (gray) are the Glauber result; the top curves (black) are the CT re- 
sult. The value of the elastic p-nucleon t-slope parameter b used in the calculation is 
indicated for each graph; VMD corresponds to b 7 v = b. 
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Figure 3.10: Integrated transparency T for A = 40, t = —2 GeV 2 , l c = 5 fm. 
The bottom curves (gray) are the Glauber result; the top curves (black) are the CT 
result. The value of the elastic p-nucleon t-slope parameter b used in the calculation 
is indicated for each graph; VMD corresponds to b y v = b. 



55 



T 

0.50 ■ 
0.45- 
0.40- 
0.35- 
0.30- 
0.25- 



T 

0.30- 
0.25- 
0.20- 
0.15- 



2 4 6 8 10 12 



(a) A = 12, g 2 = 0.5 GcV 2 , l c = 5 fm 



-t (GeV 2 ) 



2 4 6 8 10 12 



(b) A = 40, g 2 = 0.5 GeV 2 , l c = 5 fm 



-t (GeV 2 ) 



Figure 3.11: Integrated transparency T for fixed Q 2 and l c and varying t. The bottom 
curves (gray) are the Glauber result; the top curves (black) are the CT result. 
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function of Q 2 , for two different values of the coherence length. The same values of b 
and 6 7 y were used as for the T(p m = 0) calculation; VMD corresponds to b = b^y. 

The same overall features of the graphs are present as were seen for the p m = 
transparency. In addition, here one can see that for a given A and Q 2 , the trans- 
parency increases as the coherence length l c decreases, which agrees with expecta- 
tions. For the whole range of Q 2 from 2 to 12 GeV 2 , the difference between the CT 
transparency and the Glauber transparency is significant. For the higher values of 
Q 2 , the CT value is of the order of 1.5 times as large as the Glauber transparency, for 
A = 12, and 2 times as large as the Glauber transparency for A = 40. The integrated 
transparency is significantly smaller than the values for p m = 0. This is a relevant 
feature for experimentalists to note. 



In Fig. 3.11, the transparency is shown for fixed Q as a function of t. In that 



figure, Q 2 = 0.5 GeV 2 , which is small enough that for the rescattering terms (Eqs. 



3.4.7 and 3.4.8) the produced qq (at either Z2 or zs) is a normal p-meson. Thus 
no Color Transparency effects occur as it propagates from the point where it was 
produced to the point where it undergoes the hard scatter of momentum transfer 
q which knocks out the nucleon. But the large-momentum transfer scattering at Z\ 
causes the outgoing p-like configuration to be in a small-sized configuration. Hence 
the outgoing p experiences reduced interactions on its way out of the nucleus (the 
knocked-out proton also experiences reduced interactions). This is a manifestation of 
Color Transparency effects for small Q 2 (but large t). The difference between the CT 
result and the Glauber result is not as significant, however, as in the case of large Q 2 . 
The range of t shown is such that production angle of the outgoing p is small, which 
is necessary for the validity of the Glauber model. For the same Q 2 = 0.5 GeV 2 but 
for l c — 2 fm, the maximum allowable \t\ such that the p production angle is small 
is only around \t\ — 2 GeV 2 ; hence plots for this value of l c are not shown since the 
range of t would be small. 
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3.5 Conclusion 

We have calculated the transparency for 7* + A — > p + p + (A — 1)*, both without 
inclusion of CT effects (Glauber case) and with inclusion of CT effects, for several 
different combinations of A and l c . The transparencies clearly exhibit the coherence 
length effect, i.e. the decrease of the transparency as l c is increased, which is not due 
to Color Transparency. Thus to observe the effects of CT it is necessary to keep l c 
fixed while varying v and Q 2 . The quantity of experimental interest, namely the inte- 
grated transparency, is smaller in general than the transparency evaluated at missing 
momentum p m = 0. However, the difference between the Glauber transparency and 
the CT transparency is marked, particularly as Q 2 is increased while t is fixed. How- 
ever, it should still be possible to observe the effects of CT when Q 2 is small, if t is 
large enough. This represents an as yet unexplored kinematic region in the search for 
CT effects in electroproduction of vector mesons, namely small Q 2 but large t. The 
difference between the CT prediction and the Glauber prediction for the transparency 
in this case is not as large as it is in the case of large Q 2 . 
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Chapter 4 

LOW- AND INTERMEDIATE-ENERGY J/ip 
ELECTROPRODUCTION ON THE DEUTERON AT JLAB 

4.1 Introduction 

With the impending 12 GeV upgrade at JLab, electroproduction of the J/ip will be 
possible. With the mass of the J/if> being 3.097 GeV, the threshold photon energy for 
photoproduction on a single nucleon is 8.2 GeV, and is thus accessible with a 12 GeV 
electron beam. Most of the existing data on J/ip photo- and electroproduction is 
at much higher energy. The 12 GeV upgrade provides the opportunity to measure 
Jfij) production near threshold [31]. In addition, measuring electroproduction on 
the deuteron provides the opportunity to measure the J / ^-nucleon elastic scattering 
amplitude at lower energies than it has previously been measured at, if the rescattering 
of the produced J/ip on the spectator nucleon in the deuteron is significant. 

The motivation for the work in the first part of this chapter (Sees. 4.2 - 4.4) was 
a proposal at JLab [30J to measure the J / ?/>-nucleon scattering length by the reaction 
7* + d — > J/ip + p + n, where the J/ip is produced on one nucleon in the deuteron 
and then re-scatters from the other nucleon. The reason the J/^-nucleon scattering 
length is of interest is that several authors have argued that a nuclear bound state of 
the J/ip may exist [321 EH]- They propose that the force between a J/tp and a nucleon 
is purely gluonic in nature, and therefore is the analogue in QCD of the van der Waals 
force in electrodynamics, since the hadrons are of course color neutral objects. There 
is very little experimental data on elastic J / ^-nucleon scattering. There has only been 
one experimental measurement of it, at SLAC in 1977, where the J/^-nucleon total 
cross-section was extracted by measuring production of J/^'s on nuclei and using an 



59 



optical model for the re-scattering of the J/ip on the spectator nucleons [60J. 

Measurement of the scattering length provides information on the bound states of 
the two particles involved in the scattering. In particular, for an attractive potential, 
if the scattering length is positive then there exists a bound state. Since the scattering 
length is the (negative of) the zero-energy scattering amplitude, in order to measure 
this it is necessary for the two particles to scatter with small relative-momentum. In 
the case of 7* + d — > J/if> + P + n at the energies which are kinematically allowed in 
the proposed JLab experiment, it isn't possible to have an on-mass-shell nucleon and 
J /if) scatter at small relative momentum. For an incident virtual photon of energy 
v = 9 GeV, and an outgoing J/^-neutron pair with zero relative momentum, the 
minimum possible momentum of the neutron in the LAB frame (deuteron at rest) 
is ~ 0.85 GeV; for v = 6.5 GeV and zero relative momentum of the J/V'-neutron 
pair, the minimum LAB momentum of the neutron is ~ 1 GeV (see Fig. 4.4[ ). For 
zero relative momentum of the outgoing pair, the initial LAB momentum of the 
neutron (before the collision with the J/if>) must equal the final LAB momentum of 
the neutron. Therefore, the momentum of the neutron inside the deuteron would 
have to be 0.85 GeV (for v — 9 GeV). However, the deuteron wavefunction at that 
momentum is very small (essentially zero). 

So although the proposed experiment [30J may not be able to measure the J/ip- 
nucleon scattering length, it might still be possible to measure the on-mass-shell J/if>- 
nucleon scattering amplitude, but at higher relative energies. The relative energy of 
the J/^-neutron pair would still be significantly smaller than in the only existing 
data (from the 1977 experiment at SLAC). Under certain kinematic conditions, the 
dominant contributions to the amplitude will come from p-n rescattering and/or J/ip— 
n rescattering after the J /if) is produced. If we fix the magnitude of the outgoing 
neutron's momentum at a moderately large value (here taken to be 0.5 GeV) the 
contribution of the impulse diagram (where the J/if> is produced on the proton and 
the neutron recoils freely) will be negligible, since the impulse diagram is proportional 
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to the value of the deuteron wavefunction at that momentum (see Fig. 4.8 for the 
impulse and rescattering diagrams). This higher-energy rescattering is the subject of 
the second part of this chapter (Sec. 4.5). 

Note that when the relative energy of the produced particle and nucleon is small, 
the Glauber theory that was used in the previous two chapters is not applicable. In 
the Glauber theory, the projectile or produced particle always moves at high speed 
relative to the nucleons. To determine the scattering length, we have the opposite 
situation: the produced J /if) needs to be moving slowly relative to the nucleons. Also, 
in Glauber theory no account is made for the Fermi motion of the nucleons. But for 
the case of near-threshold production, the Fermi motion has a large effect on the 
amplitude and must be taken into account. Therefore a different method must be 
used to calculate the scattering amplitude in this case. For the calculations in this 
chapter, a covariant Feynman diagram method is used. 



This chapter is organized as follows. In Sec. 4J2 the diagrammatic approach is 
discussed in a heuristic manner, as well as its reduction to the Glauber theory under 



certain kinematic conditions. In Sec. 4.3, electroproduction of a particle from a nu- 
cleus is discussed. In Sec. 14.41 the kinematics for the case of zero and small relative 



momentum of the outgoing J /-?/>-neutron pair is discussed. In Sec. 4.5 the calculation 
of the invariant amplitudes for 7* + D — > J/ip+p + n are presented, including the one- 
loop diagrams corresponding to the p — n and J/^-nucleon rescattering processes. In 
order to calculate the amplitude corresponding to the low-energy J/^-neutron scatter- 
ing, which involves the scattering length, model J/^-neutron scattering wavefunctions 
and potentials are used, and it is shown that the resulting amplitude is insensitive to 
the model used. In addition, it is shown that for the kinematic conditions of the JLab 
experiment, the dominant amplitude is the impulse diagram, corresponding to J/if> 
production on the neutron with the proton recoiling freely, with no rescattering of 
any particles. This demonstrates that the measurement of the J / ^-nucleon scattering 



length is not feasible for the JLab experiment. Finally, in Sec. 4.6 calculations of the 
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amplitude for 7* + -D — > J/ip+p+n are presented under different kinematic conditions 
(not restricting the outgoing J/^-neutron pair to small relative momentum). There 
it is shown that if the J/V'-neutron elastic scattering amplitude is somewhat larger 
than the value measured at SLAC at higher energy, it may be possible to extract this 
amplitude from the Jlab experiment. 

4.2 Diagrammatic approach: Heuristic discussion 

Many years ago it was shown that the Glauber approximation can be derived using a 
diagrammatic approach ESI EZ1 ESI ES] where each Feynman diagram represents 
a Lorentz invariant amplitude, under the kinematic conditions for which the Glauber 
theory is valid, i.e. neglect of longitudinal momentum transfers, neglect of energy 
transfer from the projectile, neglect of the Fermi motion of the nucleons in the nucleus, 
and neglect of off-mass-shell effects. In this section, I'll heuristically outline this 
derivation. Here I omit various factors from expressions, in order to just show the 
dependence of quantities on the amplitudes / and wavefunction \1>. 

For example, for the case of elastic hadron-deuteron scattering, i.e. h + d — >■ 
h + d, which is the case that is treated in [35], ESI EH EH] , the Glauber result for the 
scattering amplitude consists of two types of terms, a single-scattering term and a 
double scattering term. The single-scattering term is just the product of the hadron- 
nucleon scattering amplitude and the deuteron form factor: 



where q is the 3-momentum transfer from the incident hadron h, f is the 2-body 
hadron-nucleon scattering amplitude, and G is the deuteron form factor: 



where \I/ is the deuteron momentum-space wavefunction. The second term in the 
Glauber approximation is the double scattering term, where the projectile first scat- 
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ters elastically from the proton and then scatters elastically from the neutron (or vice 
versa), and has the form 

F 2 ~ J d 2 p G(p 2 )f(q/2 + p)/(q/2 - p). (4.2.3) 

The two amplitudes F\ and F 2 can be represented diagrammatically by the diagrams 



in Fig. |4.1| If we interpret these diagrams as covariant Feynman diagrams, with the 
solid gray dot representing the deuteron vertex function T(n) (the invariant amplitude 
for the virtual dissocation D — > p + n) and the circle-with-cross representing the 
invariant amplitude for hadron-nucleon elastic scattering Ai, then the first diagram 
would have the following Lorentz invariant expression (leaving out various constants): 

d'n-, r( P )rfe)^(p,P3P,) ; (4JM) 



[n 



m 2 + ie)(p 2 — m 2 + ze)(p| — m 2 + ie 



where the factors in the denominators are due to the propagators of the internal lines, 
and n, p, etc. are 4-momenta. The second diagram would have the Lorentz invariant 
form 



,4 ,4 T(n)T(p 2 )M(p,p 3 ,k)M(n 1 fc,p 4 ) 



d n d p 2 



n 2 — m 2 + ie)(pl — m 2 + ie) , . 

i (4.2.5) 

x 

(p 2 — m 2 + ie){n 2 — m 2 + ie)(k 2 — m 2 h + ie) 

The two amplitudes F a and F& can be shown to be equal to the Glauber amplitudes 
F\ and F 2 , under the kinematic conditions where the Glauber approximation is valid: 
high-energy of the incident hadron, neglect of energy loss by the incident hadron, 
neglect of longitudinal momentum transfer at each of the /i-nucleon scatterings, and 
the non-relativistic limit |n|, | p>2 1 ^ m - hi the non-relativistic limit we have (leaving 
out various constants): 

T(p) 

2 t M ■ ~ m (4-2.6) 

and 

T(vo) . . . 
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(a) single-scattering term (b) double-scattering term 

Figure 4.1: Single- and double-scattering diagrams for hadron-deuteron elastic scat- 
tering. The internal lines with bars on them are on-mass-shell in the Glauber approx- 
imation. 



The integrals over 4-momenta are reduced to integrals over 3-momenta by integrating 
over n° and p\ using the residue theorem and keeping only the contributions due to 
the positive-energy poles at n° = \J m 2 + n 2 — ie and p\ = \/m 2 + pf— ie coming from 
the propagator denominators n 2 — m 2 + ie and p\ — m 2 + ie, respectively. Contributions 
from all other poles in n° and p\ are neglected. In F a , if the dependence of M. on the 
Fermi momentum n is neglected, one obtains 

F a ^ M J d 3 n^(n)^(n-q). (4.2.8) 

which is equal to F\ when the various constants are included. In Fb, in addition one 
more momentum component is integrated over by using 

' iir5(k 2 - m 2 h ) + V - „ 1 2 (4.2.9) 



k 2 — m\ + ie ' " k 2 — m? h 

( with V indicating the principal value) and only keeping the contribution from the 
delta function. This leaves Ff, in the schematic form 

F b ~ J d 3 nd 2 q'y*tyMM (4.2.10) 

and it is shown in (35J EE1 ETJ [38] that in the Glauber limit, F& = F 2 . 
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The Glauber approximation is thus equivalent to the sum of the covariant Feynman 
diagrams, but keeping only the positive-energy poles from some of the propagators; 
contributions from all other poles are neglected. Note that the integrations that were 
done by taking poles from the indicated propagators enforced that those particles 



are on-mass-shell in the final result. In Fig. 4.1 , those particles are indicated with 
a bar on the line. Since at each deuteron vertex one of the nucleon lines coming 
out of it is on-mass-shell, that means that the other line coming out of the deuteron 
vertex must necessarily be off-mass- shell. However, under the conditions of validity 
of the Glauber approximation, they are not very far off-shell, and the amplitudes 
Ai can be approximated by their on-shell expressions. In addition, by keeping only 
the indicated poles, the covariant diagrams become time-ordered diagrams, with a 
definite time-ordering to the sequence of hadron-nucleon scattering events. 

In the work of [351 [36l ETJ EE] the covariant expressions are further used to calcu- 
late the scattering amplitude under general conditions, when the Glauber conditions 
on momentum- and energy-transfer are not satisfied. This generalizes the Glauber 
approximation to arbitrary momentum and energy transfers. In particular it is valid 
for the case where the projectile or produced particle is moving slowly relative to 
the nucleons, which is the case we are interested in for determination of the J/ip- 
nucleon scattering length. One benefit of the diagrammatic approach with invariant 
amplitudes is that it allows one to account for the Fermi motion of the nucleons in 
the deuteron, by including the exact dependence of the scattering amplitudes M. on 
the 4-momenta of the particles; this dependence on Fermi momentum is completely 
neglected in the Glauber approximation (it assumes that the nucleons are at rest, and 
uses scattering amplitudes /(q) assuming at-rest nucleons). This becomes important 
when the scattering or particle production is near threshold, where the 2-body scat- 
tering amplitude A4 has strong dependence on the momentum of the struck nucleon. 

In [51] and [ID], the diagrammatic method is shown to reduce to the Glauber 
result for nuclei with arbitrary nucleon number A, under the Glauber conditions. In 



65 




(a) (b) 

Figure 4.2: a) Diagram included in Glauber approximation; b) Diagram not included 
in Glauber approximation. The solid black circle represents the proton-neutron scat- 
tering amplitude; the solid gray circle is the deuteron vertex function. 



addition, it is shown that the dominant diagrams that contribute to the total ampli- 
tude, under the Glauber conditions, are only those diagrams for which the nucleons 
do not interact with each other while the projectile is traversing the nucleus; diagrams 
where, e.g. 2 nucleons scatter with each other in between 2 projectile-nucleon scat- 
terings, are suppressed in the Glauber limit, compared to the other diagrams. This 
limits the number of projectile-nucleon scatterings in any diagram to at most A, and 
so there finite set of diagrams to evaluate. 

For example, for the case of incoherent scattering of a hadron by the deuteron, 
wherein the deuteron breaks up into an unbound proton and neutron in the final 



state, Fig. 4.2(a) shows a diagram which contributes in the Glauber limit (since 
the proton and neutron interact after the last projectile-nucleon interaction), while 



Fig. 4.2(b) shows a diagram which would be negligible in the Glauber limit. If the 
kinematic conditions do not satisfy the Glauber conditions, then in general both of 
these diagrams must be evaluated. 

An improvement on the Glauber approximation, called the Generalized Eikonal 
Approximation (GEA), was derived by Frankfurt, Strikman, Sargsian, et al [4T1 142| 
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EH HH SSI HE] by evaluating Feynman diagrams in the high-energy limit but retaining 
the dependence on the longitudinal momentum transfers (which is neglected in the 
Glauber approximation) and the dependence of the 2-body amplitudes on the internal 
momenta. The advantage of the GEA over the Glauber approximation is that it 
accounts for the exact momentum transfers, and therefore extends the range of validity 
beyond the Glauber limit. In addition, by virtue of its diagrammatic derivation 
the scattering amplitudes that appear at vertices have the correct dependence on 
the Fermi momentum of the nucleons. Finally, because it still assumes high-energy 
scattering, the number of Feynman diagrams that have to be evaluated is limited to 
those for which the nucleons do not interact with each other while the projectile is 
traversing the nucleus; any other diagrams are suppressed, in the high-energy limit. 
Again, this yields a finite set of diagrams to evaluate. 

The GEA would not be valid, however, for the determination of the J/^-nucleon 
scattering length, since the J/?/>- nucleon rescattering takes place at low energy, not 
high energy. Therefore in principle there would be a large number of diagrams to 
evaluate. Here we will content ourselves with only evaluating the first few diagrams. 
Previous calculations using the diagrammatic method, and keeping a finite number 
of diagrams, have been performed for 7Z) — > NNn [47] . n + D — > 7T + pn [1EIJ, 7-D — > 
n°D [49J, n~D nun [50J, np -> rjd [51] , e + D epn, e + 3 He -> epD, e + 3 He -> 
epijpn) [52], e + 3 He — > eppn [14] . — > <pD [25], and e + D — > epn [26]; the last four 
references calculate within the GEA. 

For the calculations in this chapter, it will be assumed that the diagrammatic 
expansion gives the scattering amplitude accurately, upon retaining only the first few 
diagrams (up through one- loop order, which represent the processes where the J / ip is 
produced on one nucleon and then the other (spectator) nucleon scatters afterward, 
either from the first nucleon or from the J/ip). 
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4.3 Electroproduction from a Nucleus 

We consider here electron scattering from the deuteron with production of a vector 
meson, with the final state of the proton-neutron system being a continuum state. 
The formalism for the cross-section for electroproduction from a nucleus can be found 
in Here we summarize the relevant facts. We consider here the completely 

unpolarized electron (initial and final) cross-section. Then the cross-section can be 
written in terms of the amplitudes for 7* +d — > p+n + V, i.e. vector meson production 
from virtual photons. In the LAB frame, with e' the final electron energy , Q' the 
final electron solid angle, Qy the vector meson solid angle , and p* n the final proton- 
neutron relative momentum in the p — n center-of-mass frame, the 8-fold differential 
cross-section has the form 
"V 

fi h 



de'dVL'dVL v d?p* = ( kinematic f actors ) x (v T R fi + v TT R fi + v L R fi + v TL R 

(4.3.1 



where the factors 

^ = l(/l^i(q)K)| 2 + l(/k-i(q)K)| 2 (4.3.2) 

Rf = 2Re(/iJ; i (q)ii)(/|J_ 1 (q)|i) (4.3.3) 

Rf = -2Re(/|p*(q)|i)((/|J +1 (q)|z) - (/| J_ x (q)|<» (4.3.4) 

Rf l = \(f\p(<l)\i}\ 2 , (4-3.5) 

are in terms of the matrix elements of the spherical vector components of the electro- 
magnetic current operator J between the initial deuteron state \i) and final hadron 
(p + n + J/ip) state |/). vt, vtt, etc., are kinematic factors that only depend 
on the electron momenta. R^ is the sum of the squares of the amplitudes for 
7* + d— > p + n + V for transversely polarized virtual photons, whereas R^f is an in- 
terference term between these two amplitudes. R^ is the square of the amplitude for 
a longitudinally polarized virtual photon, while R T ^ is an interference term between 
the amplitudes for production from transverse and longitudinally polarized photons. 
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The matrix element R is proportional to cos 0, while Rjf is proportional to cos 20, 
where is the angle between the plane including the initial and final electron mo- 
menta, and the plane including the 3-momentum transfer q and the J/ip momentum 
Pv- Thus if we integrate the cross-section over 0, the terms R TT and R TL drop out. 
Or, if we assume helicity conservation (i.e. the helicity of the outgoing J /if) is equal to 
the helicity of the photon) then R TT = R TL = 0. Moreover, several theoretical mod- 
els [55l [56] indicate that for small Q 2 , the amplitude for J/ijj electroproduction from 
transverse virtual photons is much larger than the amplitude for production from lon- 
gitudinally polarized virtual photons; for Q 2 = (photoproduction) the production 
amplitude for longitudinal photon polarization is of course exactly zero. Therefore in 
what follows we will neglect R^, and so the differential cross-section is simply given 
by multiplied by kinematic factors. Thus our task is to calculate 

R% = \{f\J + Mm 2 + l(/l^i(q)K)| 2 = \F + \ 2 + |F_| 2 (4.3.6) 

where F± are the amplitudes for J/ip production from positive and negative helicity 
virtual photons. In the following we will calculate the amplitude for 7* +d — > p+n+V 
by evaluating Feynman diagrams corresponding to the various processes contributing 
to it. 

4.4 Kinematics for small relative energy of the J/ip-neutron pair 

Since the scattering length is the zero-energy limit of the scattering amplitude, in 
order to measure it the relative momentum of the J/ip — n system must be small. An 
estimate of how small can be obtained by requiring only S-wave scattering, meaning 
the contribution of higher partial waves should be negligible. The classical relation 
between impact parameter and angular momentum yields an estimate for the max- 
imum I that contributes. If the relative momentum is p* and the impact parameter 
is b, then the orbital angular momentum is L = p*b = by / 2fiT* where T* is the total 
kinetic energy of the J/ip-n pair in their cm. frame, and /z is the reduced mass. The 
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largest angular momentum wave which is scattered is obtained by setting b equal to 
the range of the potential. With L 2 = 1(1 + 1) (we take H — 1), the condition for only 
S-wave scattering is that I <C 1, which implies L 2 = 6 2 2/iT* <C 1. Taking the range 
of the interaction to be ~ 1 fm yields T* <C 30MeV. 

Experimentally, perhaps the simplest quantity to measure is the total production 
cross-section, integrated over all available phase space, for a given incident photon 
energy. However, if we restrict the photon energy such that the maximum J/ip- 
neutron cm. kinetic energy T^ ax is small in of all the available phase space, then that 
means that the maximum proton-neutron relative energy will also be small everywhere 
in the available phase space: for a given value of the Mandlestam variable s for a 
system consisting of 3 particles, the total kinetic energy of any two of the particles 
(say 1 and 2) in their cm. frame satisfies 



The low-energy J /if) — n scattering amplitude is expected to be much smaller than 
the low-energy p — n scattering amplitude, and therefore the p — n rescattering would 
dominate over the J/ip — n rescattering, as contributions to the total production cross- 
section. Thus we need to restrict our considerations to a kinematic range where the 
p — n rescattering is at relatively high energy, while the J/ip — n rescattering is at 
very low energy, in order to have the possibility that the J/ip — n rescattering makes 
a noticeable contribution to the differential cross-section. 

The ideal situation would be to have the final J ftp and neutron sitting at rest in the 
LAB, with the proton moving off at high velocity. Such a final state is kinematically 
allowed for other reactions, e.g. ir + d — > npp, 7*d — > npn, but it is not possible for the 
reaction 7*rf — >■ J /if) pn, for any real or virtual photon 4-momentum. 



T12 < V s ~ m i — m 2 — m 3- 
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44-1 The T* = case 

The kinematics for the case of zero relative momentum of the J /if) — n pair is the 
simplest to analyze, since for this case the two particles are each moving with the same 
velocity (in any reference frame), and so kinematically they are identical to a single 
particle of mass M = m + my. If their common velocity in a given frame is 0, then 
the energy in that frame of the J /if) (neutron) is (3my (/3m), and the momentum of 
the J /if) (neutron) in that frame is 7/3my (7/3m). Thus the energies and momenta of 
the two particles are in the ratio Ey /E n = py/p n = ray /m. Defining Ey n = Ey + E n 
and pv n = py + p n (I use the subscript capital V to stand for the J /if) throughout 
this chapter) we have: 



TP rny 771 

— ~W^Vn 

TP m tp 

- C/ n yi - C/ Vn 

PV = IvtPVn 

Pn — mPVu 



(4.4.2) 



and also 

Ey n ~ Pvn = M 2 , (4.4.3) 
and conservation of energy and momentum is 



(4.4.4) 



v + Ed = E p + Ey n 

q + PD = Pp + PVn 

where the photon's 4-momentum is q = (u, q) and the deuteron's is Pd = (Ed, Pd)- 
Thus, just as in 2 — 2 scattering, there is only one independent variable in the final 



state, which can be taken as 9 cm (see Figure 4.3). Then we can plot the momenta of 
the proton, neutron, and J /if) in the LAB frame as a function of 9 cm . Figure 4.4 shows 
these plots for two different photon energies, one for v — 9 GeV (which is above the 
threshold energy for J /if) production on a single nucleon) and one for v = 6.5 GeV 
(which is below threshold for production on a single nucleon). 
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Figure 4.3: The reaction for T* = in the LAB and CM frames 
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Figure 4.4: Kinematics for T* = 



44.2 The T* ^ case 



Fig. |4.5| shows the minimum possible outgoing neutron LAB momentum and the 



corresponding proton LAB momentum vs. 6 cm for fixed Ty n = 30 MeV (see Fig. 4.6 
for the definition of 8 cm for the case Ty n 7^ 0); graphs are shown for photon LAB 
energy v — 9 GeV and for v = 6.5 GeV. One can see that the neutron's momentum 
is always greater than at least 0.6 GeV. Since the maximum nucleon momentum 
in the deuteron is around 0.3 GeV (the deuteron momentum-space wavefunction is 
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negligible for momenta larger than that) that means that in order for these final-state 
kinematics to occur, the neutron must have acquired its large momentum through a 
scattering event. In fact it will turn out that the dominant process corresponds to the 
impulse approximation wherein the J/ip is produced on the neutron itself, and the 
proton simply recoils freely. For the kinematics of interest here, rescattering processes 
(e.g. J/ip- neutron rescattering, J/^-proton rescattering, proton-neutron rescattering) 
make very small contributions to the total amplitude. 

We use the following notation throughout this chapter: q = (u, q) is the virtual 
photon 4-momentum in the LAB, with q 2 = —Q 2 < 0; p p is the outgoing proton 
LAB 3-momentum; p n is the outgoing neutron LAB 3-momentum; py is the J/ip 
LAB 3-momentum; and the same variables with cm superscripts denote their values 
in the overall (3-body) center-of-mass frame. 9 P , 8 n , and 9 V denote the angle that the 
outgoing proton, neutron, and J/ip momenta, respectively, make with q, in the LAB 
frame. 
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Figure 4.5: Minimum possible neutron LAB momentum (solid curve), and the cor- 
responding proton LAB momentum (dashed curve), vs. 9 crn , for Ty n = 30 MeV and 
two values of photon LAB energy v. 
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Figure 4.6: Momenta and angles in the overall cm. frame, for coplanar kinematics. 




(c) F 3a : J/ijj-n rescattering diagram 
Figure 4.7: Feynman diagrams for 7* + d — > J/ip+p + n, for production on the proton. 
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Figure 4.8: Feynman diagrams for 7* + d — > J/ip + p + n, for production on the 
neutron. 



75 



4.5 Invariant Scattering Amplitudes 



The diagrams considered here are shown in Figs. 4.7 and 4.8 There are 3 diagrams 
for production on the proton, and 3 similar diagrams where the J/ip is produced on 
the neutron. In all cases we are interested in kinematics where the J/if> and neutron 
have small relative momentum. The diagrams are covariant, and hence give Lorentz 
invariant amplitudes. In the diagrams, JvV v is the Lorentz invariant amplitude for 
the quasi-2-body process 7* + N — > V + iV (where N is a nucleon, and V stands 
for the J /if)), while Ai Vn is the Lorentz invariant amplitude for the elastic scattering 
process V + n — > V + n (with n meaning neutron), and Ai Vp and Ai np are the same 
for elastic J /if)- proton scattering and neutron-proton scattering, respectively. 



4-5.1 Impulse Diagrams 

Amplitudes F\ a and Fu are the impulse diagrams, where the J/if> is produced on 
one of the nucleons and the other nucleon (the "spectator") recoils freely without 
interacting with the other particles. In F\ a the vector meson is produced on the 
proton and the neutron is the spectator, while in Fu, the production occurs on the 
neutron and the proton is the spectator. The invariant amplitudes in this case are 

F la = M^ v (s la ,t la ) (4- 5 - 1 ) 

D(p) 

F lb = M^ v (s lb ,t lb ) T ^-. (4.5.2) 
D{n) 

Here M.^ is the Lorentz invariant amplitude for the quasi- 2-body process 7* + iV — > 
V + N (where iV is a nucleon), is the covariant vertex function for the virtual 
dissociation D — > p + n, and D(p) is the propagator denominator for the intermediate- 
state nucleon, D(p) = — p 2 + m 2 — ie. Evaluated in the LAB frame, and neglecting any 
contributions to the deuteron vertex from antinucleons, the deuteron vertex function 
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is related to the nonrelativistic deuteron wavefunction by 

^ D (k rel ) = , ~ Fd(p ) (4.5.3) 

where in the LAB frame, k re ; = p = — p n for Fi a , and k re ; = p p = — n for Fib (k re ; is 
the proton's momentum inside the deuteron, in the LAB frame, for both). 

Normalization of Amplitudes 

All 2-body amplitudes Ai lV , Ai are related to the corresponding 2-body differential 
cross-section by 

dt = 16n\(s,m 2 ,m 2 J M ^ (4 ' 5 ' 4) 

where the flux factor A is given in terms of the incident particle masses mi and m 2 
by 

A(s, m 2 , m 2 ) — (s — ml — m\) 2 — Am\m\, (4.5.5) 

and s and t are the Mandelstam variables for the 2-body process. For 7*+A^ — >■ V+N, 
we have m\ = q 2 = —Q 2 and m\ = p 2 (for F la ) or m\ = n 2 (for Fu,). Note that 
the internal nucleon lines are not on- mass-shell. For example, p 2 = {pd — Pn) 2 = 
+ m 2 — 2Mj^E n = 5m 2 — 4mE n , and so only for p„ = do we have p 2 = m 2 . 
The Mandelstam invariants that the elementary amplitudes Ai lV depend on are 
given by: 

Sla = (q + P) 2 = (Pp + Pvf 

*" = ( "~ VVY (4.5.6) 
s\b = (q + n) 2 = (p n + p v ) 2 

tu = (q-pv) 2 

Note that for the particular case of zero J/^-neutron relative momentum (T* = 0) 
we have Sx& = (p n + Pv) 2 — (m + my) 2 = M 2 where M = m + m v , so that the J/tp 
production amplitude M." iV in Fib is always at threshold. 
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In terms of the deuteron wavefunction, the amplitudes are thus: 

F la = -M^ v (s la ,t la ) ^(- Pn )v/2m(27r)3 (4.5.7) 

F lb = -M^ v (s lb ,t lb ) ii D (p p )y/2m(2nf. (4.5.8) 

The amplitudes M? v used in calculations are to be taken from experimental data on 
J/ip production on a single nucleon. 

In the above expressions for the amplitudes F la and F lb , spin labels have been 
suppressed. The initial virtual photon and the deuteron are in specific spin states, 
the final hadrons are in specific spin states, and there is a sum over the spin states 
of the intermediate-state virtual nucleon. For example, the amplitude F lb , including 
spin state specification, is explicitly: 

F lb = -y/2m(2n) 3 £ M? v {m u A, m n , Ay) ^(p p ,mi,m p ), (4.5.9) 

mi 

where mi is the spin state of the intermediate-state neutron, i.e. the line with momen- 
tum n in the Feynman diagram, m n and m v are the spin states of the final neutron 
and proton, A is the photon polarization, Ay is the J/ip polarization, and M is the 
deuteron spin state. In what follows we will assume that the elementary amplitudes 
for spin-flip are negligible compared to the non-spin-flip amplitudes, and so the am- 
plitudes are diagonal in the nucleon spin, and also in the photon and J/ip spin. In 
that case we are able to calculate the spin-averaged squares of the various amplitudes 
Fi a , F 2a , etc., and the spin-averaged square of the total amplitude. We have included 
the contribution from the D-state in the deuteron wavefunction. For v = 9 GeV 
the .D-state was found to not make a significant contribution to the amplitudes, but 
for v = 6.5 GeV the D-state did contribute significantly, especially for the impulse 
diagram Fi b . The deuteron wavefunction used was the Argonne v 18 wavefunction. 
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Parameterization of the amplitudes Ai lV 

If the cross-section for J / ip production on a single nucleon is parametrized as 

% = (4-5.10) 
at 

with the parameters A\, B\ dependent on energy (in principle), then the elementary 
production amplitude jvV v is given by 

M lV = -i v / To^A(^QV^y e |Bit (4.5.11) 

where s and t are either si a , t\ a or su, tu- 

The parameters Ai and B\ that are needed for the elementary J/ip production am- 
plitude M. lV were estimated from the (scant) existing data on exclusive J/i/j produc- 
tion on a nucleon. The only available data for the incident photon energy v ~ 10 GeV 
is from a photoproduction experiment at Cornell in 1975 [58J. For v in the range 9.3 
to 10.4 GeV, they determined A x = 1.1 ± 0.17 nb/GeV 2 = (2.8 ± 0.43) x 10~ 6 GeV~ 4 
and Bi = 1.31 ± 0.19 GeV~ 2 . Those are the values used in this analysis. 

F la and F lb 



Since F\ a is proportional to \l/(p n ), and as seen in Fig. 4.5 the outgoing neutron's 
momentum is always greater than ~ 0.6 GeV, the amplitude F\ a will therefore be very 
small, since the deuteron wavefunction is negligible for those values of momentum. 
The amplitude F^, on the other hand, is proportional to ^(p p ); thus as seen in Fig. 



4.5 for 9 cm < 0.3 rad Fu should be non-negligible for v = 9 GeV since the proton 



momentum is less than 0.4 GeV over that range of 8 cm . 



4-5.2 One-loop diagrams 



The covariant expression for a general one-loop diagram (see Fig. 4.30 ) is 

d A n V(p) M^ V M 
i(2vr) 4 D{p) D{n)D(k) 



(4.5.12) 
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Figure 4.9: General one-loop diagram, n and P2 are the same particle (either neutron 
or proton). 



where p, n, k are the internal momenta indicated in the figure, M. stands for either 
Ai pn , Ai Vn , or M Vp (elastic scattering amplitude for proton-neutron, V- neutron, or 
V-proton scattering, respectively), and D(p) = p 2 — nip + ie, etc., are propagator 



denominators. Spin labels have been suppressed in Eq. 4.5.12 in particular, there 
is an implicit sum over the spin states of the intermediate state particles (the lines 
labelled p, k, and n). There are 4 diagrams total for a given set of outgoing proton, 
neutron and J /if} momenta. Taking p 2 = p n , Pi = P P (so that the internal line n is the 
neutron, and p and k are the proton) gives one diagram (p-n rescattering diagram). 
The other 3 are: P2 = p n , Pi = py, where the internal line n is the neutron, k is 
the J /if), and p and p^ are the proton (V-n rescattering diagram); p 2 = p p , pi = p n , 
where the internal line n is the proton, and p and k are the neutron (another p-n 
rescattering diagram); and p 2 = p p , Pi = Pv, where the internal line n is the proton, 
k is the J /if), and p and p^ are the neutron (V-p rescattering diagram). 

All of the one-loop diagrams can be evaluated in the same manner; see e.g. |47j . 
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[4"5] , [30] , [51] , [IS] • We first integrate over n° by identifying the poles in the integrand 
and using the residue theorem. The contribution of a given pole corresponds to a 
particular time-ordered diagram. We make the approximation of neglecting the anti- 
nucleon contribution to the deuteron wavefunction; this corresponds to only keeping 



zero in the propagator denominator D(n) = n 2 — m 2 + ie. After doing this, and 



using the relation = x /2m(27r) 3 *(|p|) = -g^ = ^2m(2vr) 3 ^(|n|) where * is 



the momentum-space deuteron wavefunction (and we work in the rest frame of the 
deuteron), we obtain 



where we've approximated u n ~ m, which is valid since the presence of the deuteron 
wavefunction implies that only small internal momenta n contribute to the integral. 
Note that in this expression, the internal nucleon line n is now on-mass-shell, since 
n° = u n . Thus in the time-ordered diagram, only the lines p and k can be off-shell; 
all the rest are on-shell. 

The above expression for the amplitude F can be separated into two terms, one 
term in which the line k is on-mass-shell and one in which k is off-mass-shell; we follow 
closely the method of [17] . This separation is useful since the elementary amplitudes 
Ai lV , M. can be determined (at least in principle) directly from experimental data on 
the relevant 2-body scattering processes only when all 4 particles involved (2 initial 
and 2 final) are on-mass-shell. This is not the case if one of the particles involved 
is off-mass-shell. But if for reasonable choices of the off- mass- shell amplitudes the 
off-mass-shell part is small compared to the on-shell part, then the off-mass-shell part 
will not play an important role. Taking the n z axis along the direction of the vector 



the positive energy pole at n° = u n — ie, where u. 



'n — 



\/m 2 + n 2 , coming from the 




(4.5.13) 




(4.5.14) 
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Figure 4.10: Coordinate system used. For the on-shell amplitude, 6 is fixed for a 
given n = |n| and pi, P2. 



P12 = Pi + P2 gives 

D(k) = k 2 —m 2 k +ie = (pi+p 2 — n) 2 — m^-Me = Si2+m 2 — m 2 ,— 2Z? 12 u; n -|-2|p 12 |ncos6'-|-ie 

(4.5.15) 

where si 2 = (pi + P2) 2 is the Mandelstam s- variable for the elastic scattering of 
particles 1 and 2, is the mass of the real particle which the line k represents, 
E\2 = Ei + E2, and 9 is as shown in Fig. 4.10| . This allows us to write F as 



1 1 r 27V r°° r m iV m 

/ deb / dnnty(n) / dcosO - . , (4.5.16) 

7o Jo i /i 2 (n) + cos^ + 2e 



^/2m(2n) 3 2|p 12 | 
where 



/ 12 (n) = gl2 + m2 "^" 2El2 ^ (4.5.17) 
2|pi 2 |n 



is independent of 6* and 0. Now using the relation 



' -iir6(x) + V- (4.5.18) 



x + ie x 



with V representing the principal value, Eq. 4.5.16 can be written in terms of its 



on-mass-shell and off-mass-shell parts: the delta function gives the part where k is 
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on-mass-shell, since it's only non-zero (with x = cos# + f\ 2 {n)) for cos 6 = —f\ 2 {n), 
which implies k 2 = mfc- Thus we have 

F = F on + poff (4.5.19) 

where 

1 1 /"27T rn+ 

F ow = -iir — d<p dnnm{n)M lV M (4.5.20) 

v /2m(27r) 3 2|p 12 | J J\ n _\ 

and 

1 1 f 2w f°° f M lV M 

F off = r / d(j) / dnn^in) V \ dcos9—— -. (4.5.21) 

v /2m(27r) 3 2|p 12 | J r J KJ J f 12 {n) + co S v ; 

In F on , the intermediate particle line k is now on- mass-shell, so k 2 = m 2 , with k° = 
E\2 — oj n . The limits of integration |n_| and n + are the solutions of 

fn(n ± ) 2 = 1, (4.5.22) 

which are 

n± = -^L|pi 2 | ± ~^E 12 (4.5.23) 

where p 2 , E\ are the momentum and energy of outgoing particle 2 in the cm. frame 
of particles 1 and 2, and particle 2 is the same particle as the internal line with 
momentum n. Thus we have the relation 

S12 = E\ 2 - p 2 12 = (^Jm 2 +p* 2 + ^m 2 +p* 2 y. (4.5.24) 

The range of n given by |n_ | < n < n + is the range of n for which it is kinematically 
possible for the line k to be on-mass-shell (given that n, pi, and p 2 are on-shell), and 
in F on the value of cos 9 is fixed at 

a t f \ s 12 + m 2 - m 2 - 2E 12 u n 
cos6> = -fi 2 (n) = : . (4.5.25) 

The amplitudes M lV ', .M are evaluated, for a given n and 0, at that value of cos6>. 
The amplitude M, is now fully on-shell, i.e. all 4 particle lines n, k, p±, and p 2 are 
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on- mass-shell. The amplitude M. lV has only one particle off-shell (p), but given that 
the magnitude of n is small (due to the deuteron wavefunction), p is almost on-shell: 
p° = Md — uj n ~ Md — m ~ m, and so p 2 = m 2 + 0( s ^). 

In the amplitude , k is never on-mass-shell; the principal value imposes this, 
since for k to be on-mass-shell, cos# must equal — fain), which never occurs in the 
principal value. Thus the amplitudes Ai lV , M. that enter into F°^ have either one 
particle off-mass-shell (for M) or two particles off-shell (for Ai lV ). It will be shown 
below that the amplitude F°H is much smaller than F on , for the kinematics of interest 
here, and for reasonable expressions for the off-shell amplitudes A4 lV , M.. 



4-5.3 General features of the one-loop diagrams 



The on-shell part of a given one-loop amplitude is dictated largely by the behavior 



of n_ as a function of 9 cm , where for a given diagram n_ is given by Eq. 4.5.23 



A necessary condition for the on-shell part of a given one-loop diagram to be non- 
negligible is that the corresponding |n_| must be small enough so that the range of 



integration in Eq. 4.5.20 includes the momenta where the deuteron wavefunction is 



significant (see Fig. 4.11). In fact, since it is n^(n) which enters into the integral 



in Eq. 4.5.20, and this quantity is fairly sharply peaked at n ~ 0.05 GeV (see Fig. 



4.11(b)), it is necessary to have n_ < 0.05 GeV in order for the on-shell part of the 



amplitude to be non- negligible. 

For the diagram of most interest, F^ a , particles 1 and 2 are the J /if) and neutron, 



respectively. Fig. 4.12(a) shows n± vs. 6* cm for Ty n = 30 MeV, for v = 9 GeV, 



and Fig. 4.12(b) shows the same for v = 6.5 GeV, for the diagram F 3a . At both of 
these photon energies, n_ is greater than ~ 0.6 GeV for all 8 cm , and so the on-shell 
amplitude will be negligible since the deuteron wavefunction is negligible for that 



momentum. Note that for the Ty n = case (zero relative momentum of the J 'ftp — n 
pair), we have p* 2 = in Eq. 4.5.23, and so n_ = n + ; thus the on-shell amplitude, 



given by Eq. 4.5.20, is exactly zero for T Vn = 0. 



84 



¥(n) (GeV)" 3/2 
40k 

30 

20 

10 



n<P(n) (GeV)" 1/2 
0.8 i 



0.2 0.4 0.6 0.8 1 



(GeV) 




(a) 



o: 2 j ; 4 &6 l o " (GeV) 

(b) 



Figure 4.11: (a) Deuteron momentum-space wavefunction (S-state only); (b) n^{n) 
vs. n 
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Figure 4.12: n± vs. 6> cm for the amplitude F 3a (J/ip- neutron rescattering), for Ty n 
30 MeV, and v = 9.0 GeV and v = 6.5 GeV. Upper curve is n+, lower curve is n_ 
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Figure 4.13: LAB frame momenta. For a given 9 cm and T Vn , the neutron LAB momen- 
tum can range between a minimum and maximum value, both of which correspond 
to 9' = in the figure. p n y = pv + p n - 



For the amplitudes F 2a , F 2 b and Fsb, the corresponding n_ graphs are shown in 



Fig. [414] for v = 9 GeV and i/ = 6.5 GeV, for T^ n = 30 MeV. Note that at a given 
value of 8 cm and Ty n , the neutron LAB momentum can range from a minimum to 
a maximum allowed value (with these two values both corresponding to p n and py 



pointing in the same direction in the LAB, with 9' = in Fig. 4.13), and the values 
of n_ for F 2a , F 2 b and F 3 b depend on the neutron LAB momentum in addition to 
9 cm and T Vn . For our calculations, we have fixed the neutron LAB momentum for a 
given 9 cm (and T Vn = 30 MeV) at its minimum value; we denote this value by p n ,min- 
One can see from these graphs that for v = 9 GeV, there are intervals of the variable 
9 cm for which the on-shell parts of F 2a , F 2 b, and F 3b should be non-negligible, since 
n_ < 0.05 GeV there (note that for the diagram F 3 b, the J/^-proton rescattering 
occurs at relatively high energy, if the J/^-neutron relative energy is small; so is 
not directly related to the J/^-nucleon scattering length). For v = 6.5 GeV, n_ is 
larger than ~ 0.4 GeV, and so these on-shell amplitudes should be small. This is born 
out by the exact calculations, where the one-loop on-shell amplitudes for v = 6.5 GeV 



N6 




(a) v = 9 GeV (b) 1/ = 6.5 GeV 

Figure 4.14: n_ vs. 9 cm for proton-neutron rescattering and J/^-proton rescattering 
amplitudes, for Ty n = 30 MeV, p n = p n ^ min . 



are in general much smaller than those for u = 9 Ge^. 
^.5.^ Parameters used in the elementary amplitudes 



For the calculation of the amplitudes, using Eq. 4.5.20 and Eq. 4.5.21 the elementary 
amplitudes M. lV and Ai are phenomenological amplitudes obtained from existing 
experimental data. In the diffractive regime, the individual 2-body differential cross- 
sections for elastic scattering are of the form: 

da pn 



dt 

and 



A pn e Bpnt (4.5.26) 



= A Vp e B ^ (4.5.27) 

where the A's and B's can depend on energy. These are related to the elementary 
amplitudes M. by 

^ = tt^^^ orl-^l 2 (4.5.28) 

dt 16vrA(s, m\,mly 1 v ; 

where A(s, m 2 , m%) = (s — m\ — m 2 ) 2 — ^m\m\ is the flux factor. Thus we have 



M = -iJl67T\(s, m 2 , m\)A e^ Bt (4.5.29) 
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relating Ai to A and B. 

The values of A and B depend on the relative momentum (or energy) of the 



rescattering pair. Table 4.1 lists the values of the momentum p of the neutron in the 
proton's rest frame (for the p — n subsystem) and the momentum p of the J/ if) in the 
proton's rest frame (for the V — p subsystem), for the case of Ty n = 0. Note that 
these quantities are independent of 6 cm (easily shown in the overall cm. frame). Also 



given in 4.1 is the total kinetic energy of the pair in the cm. frame of that pair, T* ; 



12- 



Table 4.1: Parameters used in elementary scattering amplitudes 



v (GeV) 


Subsystem 


p (GeV) 


T* 2 (GeV) 


B (GeV- 2 ) 


a tot (mb) 


A (GeV" 4 ) 


9 


n — p 


2.25 


0.64 


5.7- 6.2 


43 - 46 


260 


9 


V-p 


7.4 


1.02 


1.31 


3.5 


1.61 


6.5 


n — p 


0.86 


0.16 


6.9 


35 


160 


6.5 


V-p 


2.84 


0.247 


1.31 


3.5 


1.61 



Given the values of the momentum p, we can determine the parameters that enter 
into the elementary amplitudes A4 pn , Ai Vp . 



p-n rescattering 

For v = 9 GeV, the existing data [53] at this momentum p give B pn = 5.7 to 6.2 GeV~ 2 . 



'pa 

The value of A pn can be obtained from the total p—n cross-section by using the optical 
theorem and neglecting the real part of the scattering amplitude: 

% I" = is* = A -> (45 ' 30) 

The measured value of a to t given in the table is then used to calculate A pn . 

For v = 6.5 GeV, the existing data give B pn = 6.9 GeV~ 2 [53J and A pn = 
160 GeV~ A Wi 
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J/ip-nucleon rescattering 

There is very little data on elastic J/^-proton scattering from which to determine the 
parameters Ay p and By p that are needed for the J/^-proton rescattering amplitude 
M Vp . For the present analysis, I have assumed the validity of Vector Meson Domi- 
nance [20] for which the t-slope for elastic J / ^-nucleon scattering is equal to the t-slope 
for the process 7* + iV -> J/tp + N, and so I've taken B Vp = B lV = 1.31 ±0.19 GeV~ 2 . 
We can obtain Ay p from the total J/^-nucleon cross-section, using the optical theo- 
rem; however, there has only been one measurement of a J t lf N [HH] • In an experiment 
in 1977 at SLAC [BO] J W photoproduction was measured on beryllium and tantalum 
targets, and the total J/^-nucleon cross-section was extracted by using an optical 
model for the rescattering of the produced J/ip from the other nucleons in the nu- 
cleus. The value they obtained was o J t ^f N = 3.5 ± 0.8 mb, which gives via the optical 
theorem Ay p = 1.61 ± 0.4 GeV~ A . In that paper, however, they also note that the 
measured J/^-photoproduction cross-section together with vector meson dominance 
arguments would give a <//■?/>- nucleon total cross-section of ~ 1 mb. So we can as- 
sume the value of the J/^-nucleon total cross-section to be not very well known. In 
addition, the photon energy in the SLAC experiment was 20 GeV, and so assuming 
forward production of the J/tp then the energy of the J/ip in the LAB frame would 
also be ~ 20 GeV, giving a kinetic energy in the LAB of ~ 17 GeV. This is signifi- 
cantly larger than the kinetic energy of the J/ip in the proton rest frame considered 
here, where for v = 9 GeV it is 4.94 GeV and for v = 6.5 GeV it is 1.1 GeV. This 
introduces more uncertainty in the value of Ay p to be used. In [33], a theoretical cal- 
culation of the J/V>-nucleon scattering length yields a value for the total J/^-nucleon 
cross-section at threshold of 7 mb, and it is argued that the total cross-section should 
decrease as the energy is increased from threshold. Thus at the energy of the J/ip- 
proton rescattering here, the value of Ay p may be larger than the value measured in 
the experiment at SLAC. For the purpose of calculating the amplitude F 3 b, however, 
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we will use the value measured at SLAC. (In Sec. 4.6 on intermediate-energy J /if) 
production, the uncertainty in the value of Ay v will be put to good use, as Ay v will 
be considered a free parameter, which might in fact be measured by the proposed 
experiment at JLab). 

4-5.5 Subthreshold J /if) production 

The threshold photon energy for production on a single nucleon at rest is 

"thresh = m v H (4.5.31) 

2m 

while for production on the deuteron it is 

m 2 v + Q 2 

"thresh = m v H — — — (4.5.32) 

ZlVld 

For Q 2 = 0.5 GeV, these are 8.47 GeV and 5.78 GeV, respectively. 

For v = 6.5 GeV, which is below threshold for J /if) production on a single nucleon 
at rest, we assume that the production mechanism is the same as for production on 
a free nucleon. The Fermi motion of the nucleon in the deuteron is what allows the 
production to occur, e.g. if the nucleon is moving towards the photon with a large 
enough momentum then the value of Si = (q + p) 2 , where p is the 4-momentum of 
the nucleon in the deuteron, will be above the threshold value. In the calculation of 
the amplitudes for v = 6.5 GeV this condition was imposed on the internal nucleon 
momentum in the integrals involved. 

4-5.6 Calculated On-shell and Off-shell amplitudes 



Using the parameters in Table 4.1, the on-shell and off-shell parts of the amplitudes 



were calculated. The squares of the individual amplitudes F 2a , F 2 b, and F 3 b are shown 



in Figs. 4.15 and 4.16 shown in the graphs is a curve which includes only the (square 
of the) on-shell part of the amplitude, and also a curve which is the square of the total 
amplitude including both the on-shell and off-shell parts. For the off-shell parts, the 
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same parametrizations of the elementary amplitudes JvV v and M. were used as for 
the on-shell parts. As stated previously, the off-shell parts are very small compared 
to the on-shell parts, which means that knowledge of the exact forms of the off-shell 
elementary amplitudes M. lV and A4 are not needed. 

Since the J / -^-nucleon scattering length is expected to be small (much smaller than 
e.g. the proton-neutron scattering length), the J/^-neutron rescattering diagram F 3a 
should be a small contribution to the total amplitude. This is born out in the next 
subsection, where F 3a is calculated using a model potential and wavefunction, for a 
value of the scattering length of the order of that predicted by theoretical models. 

4-5.7 J /ip -neutron Rescattering Diagrams and the Scattering Length 



Diagram F 3a (Fig. 4.17) is the J/V'-neutron rescattering diagram. This is the diagram 
where the J/ip and neutron scatter from each other with small relative momentum; 
hence this amplitude will involve the scattering length for the J / ?/>-neutron interaction. 

The amplitude F 3a 

The on-shell and off-shell parts of F 3a are 

1 1 ftn pn+ 

F£ = -m—===- - dcf> dn n^(n)M^ v M Vn (4.5.33) 

^2m(2iry ^\Pvn\ Jo J\ n _\ 

and 



ft 1 1 r 2n f°° f M lV M Vn 

F° f J = —== , r / d(j) / dnn^in) V \ dcos9—— -. (4.5.34) 

3a v /2m(27r) 3 2|py ri ,| J Y J KJ J / 12 (n) + cos0 v ; 



where py n = Pv + Pn; an d is the angle between p Vn and n (see Fig. 4.10). For 



T* n = we have n+ = |n_| (see Eq. 4.5.23), and so F 3 ° a n = 0. For T^ n = 30 MeV 



Fg" will be small because n_ > 0.6 GeV for the possible JLab kinematics . Thus the 

main contribution to F 3a is from F£*, for which the intermediate-state J/ip is always 

off-mass-shell. However, because of the propagator denominator fain) + cos 6, where 

f i \ , a s Vn + m 2 -m 2 v -2E Vn u n 1 / 2 2 \ 

h2\n) + cos = -. hcosfc' = - — (k - m v , (4.5.35) 

2p Vn n 2pv n n\ / 
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\F 2 a\ 2 (GeV" 2 ) 




(c) 

Figure 4.15: Squares of amplitudes for v = 9.0 GeV, Ty n = 30 MeV, p n = p n , m in- 
The solid curve is the total (on-shell plus off-shell parts), while the dashed curve is 
only including the on-shell part of the amplitude. 
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Figure 4.16: Squares of amplitudes for v = 6.5 GeV, Ty n = 30 MeV, p n = p n ,min- 
The solid curve is the total (on-shell plus off-shell parts), while the dashed curve is 
only including the on-shell part of the amplitude. 



93 




D 



11 



Figure 4.17: F 3a 



contributions to F^ from values of n for which k is far off-mass-shell will be small. 
To obtain estimates of , we will therefore evaluate it using on-mass-shell values 
of M^ v and M Vn . 

For small relative momentum of the J/^-neutron pair, Ai is related to the J/ip- 
neutron scattering length. The relation between the invariant amplitude Ai Vn and 
the scattering amplitude f(k,9) is [63] 



for the on-energy-shell amplitudes. Note that in the amplitude Ai is off- 

energy-shell if it's on-mass-shell, since the magnitude of the final J/V'-neutron relative 
momentum is not equal to the magnitude of their initial relative momentum, since n 
is being integrated over. 

Scattering Length 

The definition of the scattering length for a given 2-body interaction is that it is 
the (negative of) the zero-energy limit of the scattering amplitude in the center of 
mass frame. For low-energy scattering, only the S-w&ve will contribute, and so the 
scattering amplitude is just a constant (independent of scattering angle), i.e. 



M 



87iyi^/(M) 



(4.5.36) 



lim/(fc,0) 



= —a 



(4.5.37) 
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which defines the scattering length a [65J. The scattering amplitude / here is the 
on-energy-shell amplitude, so the initial and final relative momenta are k, and k/ 
with |kj| = |kj| = k. We require the off-energy-shell amplitude, which is given by 

/ w (k!,k a ) = -(2tt)V (k 2 |V|< } ) = -(2tt)V K; } |V|ki) = -(2vr)V (k!|V|<V 

(4.5.38) 

where \l is the reduced mass; \ti is the initial relative momentum (in terms of n, k in 



Fig. 4.17); k 2 is the final relative momentum (in terms of py, p„ in Fig. 4.17); and 



^/k 2 is the exact scattering wavefunction for asymptotic relative momentum k 2 . Since 
this off-energy-shell scattering amplitude depends on the scattering wavefunction 
and the potential V, both of which are unknown for J/^-nucleon elastic scattering, 
we will resort to models in order to estimate the amplitude. 

We normalize our S"-wave wavefunction and define the radial wavefunction u(r), 

by: 

#( r ) = _^ == e m)^l (4.5.39) 

where 5(h) is the S'-wave phase shift. In order to calculate the matrix element, 
we can either specify a model potential and solve the Schrodinger equation for the 
wavefunction u(r), or instead specify a model zero-energy wavefunction u°(r) which 
determines the potential V(r) via the Schrodinger equation, and use that potential 
to solve for the wavefunction for k ^ 0. We will do the second procedure, choosing a 
model zero-energy wavefunction which satisfies the minimal constraints imposed by 
the Schrodinger equation. 

We assume the J/^-nucleon potential is of finite range, and so is zero for r larger 
than some distance R. The phase-shift S(k) satisfies the following well-known prop- 
erties |65] as k — >• : 



1. for a repulsive potential, or an attractive potential that doesn't admit a bound 
state: 5 — > —ak as k — > 0; 
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2. for an attractive potential which admits a single bound state: 5 — > n — ak as 
k ->■ 

The zero-energy wavefunction \I/(r) outside the range of the potential is then 

for both cases, while the zero-energy radial wavefunction u° out differs by a minus sign 
for the two cases; this is purely due to including the factor e i5 ^ k ' in the definition of \I> 



in Eq. (4.5.39) (see Appendix D for more details). The superscript indicates k = 0. 
Furthermore, my normalization conventions give a > for either a repulsive potential 
or an attractive potential with a bound state, and a < for an attractive potential 
that doesn't admit a bound state. In all cases a is the intercept on the r-axis of 
u° out (r) (Fig. p8b. 



Theoretical calculations [221 El] give values of \a\ ~ 0.3 fm, with effective range 
r e ~ 2.0 fm. It is thought that the potential is attractive, but too weak to support 
a bound state. Below, calculations of F^ a are made for both cases of an attractive 
potential: a > (bound state) and a < (no bound state). 

4-5.8 Positive scattering length a 

For the case of a positive scattering length and attractive potential (which possesses a 
bound state), we have u® ut (r) = — (r — a) . The properties imposed on the zero-energy 
wavefunction u for r < R, required by the Schrodinger equation, are then: 

u in {0) = (4.5.41) 

uUR) = u out (R) = -{R - a) (4.5.42) 
< n (R) = <ut(R) = -1 (4-5.43) 
< n {R) = (4.5.44) 
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(a) Attractive potential that possesses a bound state, a > 
in this case. 
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a < 




V(r) 



(b) Attractive potential with no bound states, a < in this 
case. 




a > 

(c) Repulsive potential, a > in this case. 



Figure 4.18: Scattering length, u out {r) (dashed line), and Ui n (r) (solid curve) for 
different types of potentials. 
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The first 3 are the standard requirements imposed by the Schrodinger equation, while 
the last one IS the Schrodinger equation at r = R because V(R) = 0. We will assume 
the simplest form for u that's consistent with these requirements, which is that u is 
a cubic function of r. 



Model wavefunction for k = 



The wavefunction u in that satisfies Eqs. 4.5.41 - 4.5.44 is 
u° in (r) = (-1 + -^)r + f-^SV 2 + 



RZJ 



(4.5.45) 



RJ' ' V R?) 

One further requirement on u\ n is that u® n have no zeros on the interval [0, R] (besides 
at r = 0). This ensures that the corresponding potential V(r) is non-singular, since 
for k = 0, V(r) = — . The zeros of are at r = and 

' \ ' 2/Lt u ln 



R 2 2V a 

and the right-hand-side must lie outside the range [0, 1]. This requires either 



(4.5.46) 



R< a 



(4.5.47) 



or 



R > 3a 



The potential is 



V(r) 



1 1 



r-R 



2firR 2 (3- f) -3Rr + r~ 



(4.5.48) 



(4.5.49) 



and one can see that for R > 3a the potential is repulsive. Therefore we require 
R< a. 



Model wavefunction for k ^ 

Given this model potential we can proceed to calculate the off-shell scattering am- 



plitude Eq. ( 4.5.38 ) and the amplitude F 3 ° a // , Eq. ( 4.5.34 ), once we calculate the 
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(a) u m (r) for k = 0, 50, and 100 MeV. 
Solid black curve is for k = 0, solid gray 
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(c) Off-shell amplitude for k — (solid 
black), 50 (black dashed), and 100 MeV 
(gray dashed) 
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Figure 4.19: Wavefunction and off-shell amplitudes for a = 0.3 fm, R = 0.1 fm, 



using the model potential of Eq. 4.5.49 



wavefunction u for non-zero k for a given a and R. The potential does not admit 
analytic solutions for non-zero k, so we calculate them numerically. However, we only 
require them for small k (less than 100 MeV or so). Taking a = 0.3 fm, R = 0.1 fm, 



Fig. |4.19(a) and Fig. 4.19(b) show u in for k = 0, 50 MeV, and 100 MeV; there's 



virtually no difference between them. Figs. 4.19(c) and 4.19(d) show the off-energy- 
shell amplitude f Vn (p,k) for k = 0, 50 MeV, and 100 MeV; again there's virtually 
no variation of f Vn {p, k) for k up to 100 MeV. 
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Since we are only interested in the J/^-neutron relative momentum up to around 
100 MeV, it is legitimate to approximate the off-energy-shell amplitude f Vn (p,k) ~ 
f Vn (p, 0) for the range of k we are interested in. For our model wavefunction u° in we 
can evaluate f Vn (p,0) analytically: 

r n (p,o) = -(2tt)v (pimj-y = -27T 2 (p\u\<- ] r 

— 2n 2 f 1 r R 

The momentum p appearing in f(p, 0) is the relative momentum of the J/^-neutron 
pair in their center-of-mass frame, before they scatter in diagram F 3a ; it is thus the 
magnitude of n (or -k) in the outgoing V — n center-of-mass frame, and so we must 
boost n to that frame. 

4-5.9 Negative scattering length a 

If the potential is attractive but too weak to support a bound state, then a < and 
we have u° out (r) = r — a . The properties imposed on the zero-energy wavefunction u 
for r < R, required by the Schrodinger equation, are then: 

«in(0) = (4.5.51) 
u in (R) = u out {R) = R-a (4.5.52) 

< n (R) = <ut(R) = 1 (4-5.53) 
u'i{R) = (4.5.54) 
In this case the zero-energy wavefunction is 

^(r)=(l-|)r+(|)r 2 -(A)r 3 . (4.5.55) 

The requirement that u have no zeros on [0, R] imposes no restriction on a and R 
in this case. Theoretical calculations give a around —0.3 fm, and effective range 



100 



w in (r) (fm) 
1.5 




0.2 0.4 0.6 0.8 1.0 1.2 



r(fm) 



(a) u m {r) for k = 0, 50, and 100 McV. 
Solid black curve is for k = 0, solid gray 
curve is for fc — 50 MeV; dashed curve is 
for k = 100 MeV. 
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(b) Off-shell amplitude for k = (solid 
black), 50 (black dashed), and 100 MeV 
(gray dashed) 



Figure 4.20: Wavefunction and off-shell amplitudes for a = —0.3 fm and r e = 2.0 fm 



for the model potential Eq. 4.5.49 



2.0 fm 



Using these values with our model wavefunction and potential 



implies R — 1.3 fm. The wavefunction and off-energy-shell scattering amplitude for 



this case are shown in Fig. 4.20 Again there's very little difference between the curves 
for k from to 0.1 GeV, so to calculate F 3a in this case we will again approximate 
f Vn (p, k) ~ f Vn (p, 0). We have for a < 0: 



r n (p,o) 



6a / sin pi? 



p 2 R 2 



pR 



(4.5.56) 



4-5.10 Results 



The results for all amplitudes are shown in Figs. 4.21 and 4.22 for positive a, and 
Figs. 4.23| and 4.24 for negative a, for the case of T* = 0.. As can be seen in the 
graphs, the amplitude F^ a is smaller than F^, but in the range of 9 cm where these two 
diagrams (the J/^-nucleon rescattering diagrams) dominate, the square of the total 
amplitude is much smaller than the square of the total amplitude at its maximum, by 
a factor of 10 5 to 10 3 . The total amplitude is due almost exclusively to the impulse 
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(a) Squares of individual amplitudes. Dashed curve is the square of the total 
amplitude. Not shown are F\ a , F^a, which are negligible. 
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Figure 4.21: Squares of amplitudes for v = 9 GeV, T* = 0, and a = 0.3 
R = 0.1 fm. 
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(a) Squares of individual amplitudes. Dashed curve is the square of the total 
amplitude. 
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(b) Square of total amplitude. Solid: includes all amplitudes. Dashed: omitting 



Figure 4.22: Squares of amplitudes for v = 6.5 GeV, T* = 0, and a = 0.3 
R = 0.1 fm. 
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(a) Squares of individual amplitudes. Dashed curve is the square of the total 
amplitude. Not shown are F\ a , i^a, which are negligible. 
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(b) Square of total amplitude. Solid: includes all amplitudes. Dashed: omitting 
F 3a (not distinguishable from the solid line at this scale). 



Figure 4.23: Squares of amplitudes for v = 9 GeV, T* = 0, for a = —0.3 
R = 1.3 fm, r e = 2.0 fm, for the model wavefunction of Eq. (4.5.55). 
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Figure 4.24: Square of total amplitude for v = 6.5 GeV, T* = 0, for a = —0.3 fm, 



R = 1.3 fm, r e = 2.0 fm, for the model wavefunction of Eq. (4.5.55). Solid: includes 
all amplitudes. Dashed: omitting F 3a . 



diagram F^, in the region where the total amplitude is non- negligible. 

F 3a was also calculated using 3 other potentials, for T* = 0: a square- well potential 



yielding a = 0.3 fm and R = 0.1 fm, and also the potential of Eq. 4.5.49 but with 
a = 0.3 fm, R = 0.29 fm, and a square-well potential yielding a = 0.3 /m but with 
R = 0.29 fm. The results are shown in Fig. |4.25| It can be seen that there's not 
much difference in the value of F 3a for different potentials with the same scattering 
length. For v = 6.5 GeV, there's virtually no difference at all. 



4-5.11 Non-zero T* and Differential Cross- sections 

The electroproduction differential cross-sections were calculated for the case of T* = 
30 MeV. Coplanar kinematics were used (i.e. all 3 final momenta p p , p n , and py he 
in the same plane). At a given value of the proton momentum p p , for a fixed value 



105 



f (p,k) (fm) 




TSbP (fm_1) 



(a) Off-shell amplitude for k = for a 
square-well potential, for v = 9.0 GeV, 
with a = 0.3 fm, R = 0.1 fm (the depth 
is U = 13.13 GeV- 2 = 337.3 fm 2 ). 
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( c ) |-P3a| 2 for 4 different potentials, for 
v = 6.5 GeV. Top curve is for both the 
cubic wavefunction and the square- well for 
a = 0.3 fm, R — 0.1 fm. Bottom curve is 
for both the cubic wavefunction and the 
square-well for a = 0.3 fm, R = 0.29 fm. 
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(b) |-F3a| 2 for 4 different potentials (y = 
9.0 GeV). From top to bottom the curves 
are: cubic wavefunction with a = 0.3 fm, 
R = 0.1 fm; square- well potential with 
a = 0.3 fm, R = 0.1 fm; cubic wavefunc- 
tion with a = 0.3 fm, R = 0.29 fm; square- 
well potential with a = 0.3 fm, R = 0.29 
fm. 
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Figure 4.25: Off-shell amplitudes for various potentials, for T* = 0. 
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of Ty n , the magnitude of the neutron's momentum in the LAB satisfies 

|.Pn,rran| — Pn — Pn,max (4.5.57) 

where 

E* p* n 

Pn,min = —==\pVn\ -==Ey n (4.5.58) 

and 

Pn,min = -^=\p Vn \ + ~y^E Vn (4.5.59) 
V s Vn V s Vn 

and where £"* is the neutron energy in the J/^-neutron cm. frame, p* n is the neutron 
momentum in the J/^-neutron cm. frame, Ey n is the total energy of the J/^-neutron 
pair in the LAB frame, py n is the total momentum of the J/^-neutron pair in the 
LAB frame, and sy n = (py +p n ) 2 - Note that p n ,min and p n , max depend on p p , for a 
given value of Ty n . 

For the differential cross-section calculated here, at a given value of p p the value of 
the magnitude of the neutron's momentum was set at the p n ,min corresponding to that 
Pp. This corresponds to the neutron and J/ip momenta both pointing in the same 
direction in the LAB frame (collinear momenta); in the J/ip- neutron cm. frame, the 
neutron momentum points opposite the direction that the total neutron plus J/ip 
momentum points in the LAB, while the J/ip momentum in the J/ip- neutron cm. 
frame points in the same direction that the total neutron plus J/ip momentum points 
in the LAB. 

The electroproduction differential cross-section in the LAB frame is given by 

d ± = V ° Vt E ' x ^- V l ^ Iff (4.5.60) 

dE'dVL'dppdVtpdVLn 8(2?r) 3 M d E 8(2tt) 5 E p \E v p 2 n - E n p n ■ p v \ 1 1 v 1 

where F = F la + F lb + F 2a + E 2 b + F 3a + F 3b is the total amplitude for J/ tp production 

from a virtual photon, and 



^1QE 2 E' 2 - Q 4 = 4EE'cos 2 (9'/2), (4.5.61) 
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and 



E' = E - v = 12 GeV - v. 



(4.5.63) 



In the above, E is the initial electron energy (taken to be 12 GeV), E' is the final 
electron energy, and 9' is the scattering angle of the electron relative to the initial 
electron momentum (all quantities in the LAB frame). 

The results are shown in Figs. 4.26, 4.27 and 4.28 for the case of Ty n = 30 



MeV. For v = 9 GeV, the squares of the individual amplitudes are shown in Fig. 
4.26(a)| along with the square of the total amplitude including all diagrams. As 
can be seen in that graph, the amplitude F 3a makes a negligible contribution to the 
total amplitude. There are intervals of 9 cm where the amplitudes F lb , F 2 b, and F 3b 



individually dominate the total amplitude. However, by comparing Fig. 4.26(a 



with Fig. 4.26(c), which shows the electroproduction differential cross-section on a 



linear scale, over the range of 6 cm for which the cross- section is non- negligible (for 
< cm < 0.2 rad) the cross-section is due exclusively to the impulse diagram F lb . 
(The very small "bump" visible in Fig. 4.26(c) at 9 cm ~ 0.3 rad is due to the proton- 
neutron rescattering amplitude F 2 b-) 

For v = 6.5 GeV, the difference between the cross-section including F 3a and omit- 



ting F 3a is visible in the logarithmic-scale graphs (Figs. 4.28(a) and 4.28(b)), but 



not in the linear-scale graph, Fig. 4.28(c) However, by comparing Fig. 4.27 with 



Fig. 4.28(c), one can see that over the range of 9 cm for which the cross-section is 



non- negligible (for < 9 cm < 0.6 rad) the cross-section is due exclusively to the im- 
pulse diagram Fi b . Thus the rescattering effects are negligible, for these kinematics, 
for both v = 9 GeV and v = 6.5 GeV. 



4-5.12 Conclusion 

It does not appear to be possible to measure the J/^-nucleon scattering length via 
production on the deuteron, under the kinematic conditions available at JLab. For 
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(a) Squares of individual amplitudes, for positive a. Dashed 
curve: square of total amplitude. Impulse: Fit,. P ~ n: Fzb- 
V — p: Fsb. V — n: F^ a . Not shown are F\ a , i^a, which are 
negligible 

nb/GeV 2 




0.5 1.0 1.5 2.0 2.5 3.0 



0cm (rad) 



(b) Electroproduction differential cross- 
section (logarithmic scale). 
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(c) Electroproduction differential cross- 
section (linear scale). 

Figure 4.26: Squares of amplitudes, and electroproduction differential cross-section 
Eq. (4.5.60), for v = 9 GeV, Ty n = 30 MeV, p n = p n ,min- Solid curves in (b) and 



(c) includes all amplitudes, dashed curves (not distinguishable from the solid curve) 



omit F 3a , for the model potential of Eq. (4.5.49) 
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Figure 4.27: Squares of individual amplitudes for v = 6.5 GeV, T Vn = 30 MeV, 
Pn = Pn,min- Dashed curve is the square of the total amplitude, a = 0.3 fm. 



small values of the relative momentum of the outgoing J/^-neutron pair, the ini- 
tial momentum of the neutron inside the deuteron that is required for on-mass-shell 



rescattering of the J/i/j- neutron pair is larger than ~ 0.6 GeV (see Fig. 4.12), where 



the deuteron wavefunction is negligible. The off-mass-shell part of the rescattering 
amplitude was calculated using model J/^-nucleon potentials and was found to make 
a negligible contribution to the total amplitude. The vast majority of J/ip production 
events, for Ty n < 0.03 GeV, will be at small values of 9 cm , where the impulse diagram 
Fib dominates, and therefore information on J/^-nucleon elastic scattering at small 
relative energy cannot be obtained. 



4.6 Intermediate energy J/ip production on the deuteron 



It may be possible to extract the J/^-nucleon elastic scattering amplitude from the 
7* + -D — > J ftp + p + n experiment, at higher relative energy of the J/^-nucleon pair, 
under different kinematic conditions for the final-state particles than was considered 
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Figure 4.28: Electroproduction differential cross-section Eq. (4.5.60) for v = 6.5 GeV 



Ty n = 30 MeV, p n = Pn^mini and a = ±0.3 fm. In (c), the dashed curves (not 
distinguishable from solid curve) are omitting F^ a for both positive and negative a, 



for the model potential of Eq. (4.5.49). 



Ill 



in the previous sections of this chapter. Under certain kinematic conditions, the 
dominant contributions to the amplitude will come from rescattering diagrams (p- 
n rescattering and J/ip — n rescattering). If we fix the magnitude of the outgoing 
neutron's momentum at a moderately large value (here taken to be 0.5 GeV) the 
contribution of the impulse diagram will be negligible, since the impulse diagram is 
proportional to the value of the deuteron wavefunction at that momentum (see Fig. 



4.8 for the impulse and rescattering diagrams). For the analysis presented here, we: 

• use coplanar kinematics 

• fix the magnitude of the outgoing neutron momentum at p n = 0.5 GeV 

• fix the 4- momentum-transfer-squared t = (q — pv) 2 at a particular value 

• plot amplitudes or differential cross-sections vs. 9 n (the angle that the outgoing 
neutron momentum p n makes with direction of the incoming photon momen- 
tum) for fixed p n and t 



see Fig. 4.29). For some range of v and t, these graphs will display peaks due to 



p — n and J/ip — n on-mass-shell rescattering. For v — 10 GeV, the peak due to 



J/ip — n rescattering is evident (see Fig. 4.38), but for v — 9 GeV it is not evident 



(see Fig. 4.40). This analysis is similar to what has been done in |47J for the reaction 

The kinematics here are very different than in the previous sections. There it was 
the relative energy of the J/^-neutron system which was kept fixed, at a small value, 
while the parameter which was varied was the angle of the proton momentum in the 
overall center-of-mass system. 
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Figure 4.29: General orientation of momentum vectors, for coplanar kinematics, q is 
the photon momentum vector. 



4-6.1 Intermediate- energy J/ip production 

We are interested here in kinematics available at JLab after the 12 GeV upgrade. 
The maximum (virtual) photon energy is then around 11 GeV. Here we evaluate 
the amplitude for virtual photon 4-momentum q = [y, q) with v = 10 GeV, Q 2 = 
—q 2 = 0.5 GeV, keeping p n = 0.5 GeV and t fixed. The results presented here are 
for t = —2 GeV 2 ; calculations were done for larger values of \t\, with similar results 
(although the total amplitude decreases with increasing \t\). We consider the same 



set of diagrams as before, shown in Figs. 4/7 and 4^ The impulse diagrams, F\ a and 
Fu, are negligible for these kinematics. 

One-loop diagrams - General features 



Following [47], the main features of the on-shell part of the one-loop amplitudes in 
this case can be seen by first approximating the elementary amplitudes Ai lV and A4 
to be constants, evaluated at a typical value of the intermediate-state momentum n 
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Figure 4.30: General one-loop diagram, n and p<i are the same particle (either neutron 
or proton). 



(Fig. 4.30). In that case we obtain for the on-shell amplitude: 

1 2tt 



—ITl- 



M lV M / dnn^(n) 



(4.6.1) 



|n_| 



y/2m{2ir) 3 2|pia | 

For a given value of |p n | and t, most of the dependence on 9 n in the above expression 
comes from the dependence of n_ on 6 n : if |n_| is larger than around 0.4 GeV, the 
amplitude is essentially zero since the deuteron wavefunction is essentially zero for 
n > 0.4 GeV. The maximum of F on as a function of 9 n occurs at the value of 9 n for 
which n_ = 0. The momentum |p»i2 1 also varies with 6 n , but it is much more nearly 



constant than n_ is (Figs. 4.31 and 4.32). 

Note that for a given t, p n , and 6 n , there are 2 sets of allowed values of the proton 
and J/ip momentum {p p , py}; I've called the two sets the "plus" set and the "minus" 



set. If we define x and z axes as in Fig. 4.33, with the x-component of the neutron 



momentum always positive, then the "plus" kinematics is as shown in Fig. 4.33(a) 



and the "minus" kinematics is as shown in Fig. 4.33(b) As seen from Fig. 4.34 



for the "plus" kinematics, p px is negative for all 9 n (while py x takes both positive 
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Figure 4.31: pi 2 vs. n , for photon energy v — 10 GeV, and £ = — 2 GeV 2 . (a) and 
(b) are for the "plus" kinematics, (c) and (d) are for the "minus" kinematics. 
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Figure 4.32: py p vs. 9 n , for photon energy v — 10 GeV, and £ = — 2 GeV 2 . (a) is for 
the "plus" kinematics, (c) is for the "minus" kinematics. 
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and negative values over the range of 9 n ), while for the "minus" kinematics, py x is 
negative for all 9 n (while p px takes both positive and negative values over the range 
of 0A 



Figs. 4.35 and Fig. 4.36 show graphs of |n_| vs. 9 n for the 3 different pairs of 
outgoing particles. Since the value of |n_| varies greatly with 9 n , the value of the 
integral Jj™ + | dn n^{n) varies greatly also, and has a prominent peak at the value of 



9 n for which |n_| = 0. This can be seen from Fig. 4.37, which shows the integral 
vs. 9 n . The general shape of these graphs is retained in the calculation of the exact 
amplitudes (i.e., including all dependence of the elementary amplitudes M? v and M. 
on the internal momentum n), including the location of the peaks. 

Calculation of Amplitudes 

For the calculation of the amplitudes, the elementary amplitudes Ai lV and hA are 
taken to be of the diffractive form Ae^ Bt with parameters determined from existing 
experimental data. For the <//■?/>- nucleon rescattering diagrams, the only available 



data is from the experiment at SLAC [60j discussed in Sec. 4.5.4 They determined 
the total J/^-nucleon cross-section to be a J t lf N = 3.5 ± 0.8 mb, which gives via the 
optical theorem Ay n = 1.61 ± 0.4 GeV~ A . The energy of the J/ip in this experiment 
was ~ 20 GeV in the Lab frame (nucleon at rest). However, for our kinematics 
the rescattering of the J/ip on the nucleon takes place at an energy in the outgoing 
neutron's rest frame of from 6 to 10 GeV, which is significantly smaller than in the 
SLAC experiment; thus the value of Ay n at our energy may be significantly different. 
Since the entire reason for measuring the cross-section for this process is to extract the 
J/^-nucleon scattering amplitude in an energy region where it has not been measured 
before, I've used several different values of the parameter Ay n in the calculations, from 
the value measured at SLAC up to 10 times the SLAC value. Since the total cross- 



section a tot for J/^-nucleon scattering goes like y/ A Vn (Eq. 4.5.30), this corresponds 



to a range of a to t (which is what was actually measured in the SLAC experiment) 
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Figure 4.33: (a) "plus" kinematics and (b) "minus" kinematics. For "plus", p p is 
always on the opposite side of the photon momentum q direction as the neutron 
momentum. For "minus" , py is always on the opposite side of the photon momentum 
q direction as the neutron momentum. 
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Figure 4.34: (a) p p:r vs. # n and (b) py x vs. n , for photon energy v — 10 GeV, and 
£ = — 2 GeV 2 . The solid curves are the "plus" kinematics, and the dashed curves are 
the "minus" kinematics. 
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Figure 4.35: vs. 0„, for photon energy v — 10 GeF, and t — —2 GeV 2 . (a) and 
(b) are for the "minus" kinematics, (c) and (d) are for the "plus" kinematics. 



|n_| (GeV) 




~20 40 60 W~ 6n (deg) 

(a) J/ip-p rescattering 



|n_| (GeV) 




(b) J/^-p rescattering 



Figure 4.36: vs. 9 n for J/ip-p rescattering, for photon energy v — 10 GeV, and 
t = — 2 GeV^ 2 . (a) is for the "minus" kinematics, (b) is for the "plus" kinematics. 
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Figure 4.37: The integral from Eq. (4.6.1), for photon energy v = 10 GeV^, and 
t = —2 GeV 2 . (a) and (b) are for the "plus" kinematics, (c) and (d) are for the 
"minus" kinematics. 
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of from 1 to ~ 3 times the SLAC value. In [23], theoretical calculation of the J/if>- 
nucleon elastic scattering cross-section at threshold yielded 7 mb, which is twice the 
value measured at higher energy at SLAC. 

The full calculation of the amplitudes must of course include the off-shell parts. 
If we use the same parametrization of the elementary amplitudes as in the on-shell 
part, then the off-shell parts are very small compared to the on-shell parts. Fig. 



4.38 shows \F 2a + F 3a \ 2 (which are the amplitudes for production on the proton and 
rescattering on the neutron) as a function of 9 n . The negative values of 9 n are for 
the "minus" kinematics, and the positive values are for the "plus" kinematics. It can 
be seen that the difference between the on-shell amplitude and the total (on- plus 
off-shell) amplitude is negligible. Three different values of the J/V'-nucleon amplitude 
coefficient Ay n are used. It is seen that only if Ay n is of the order of 10 times as large 
as the previously measured value is there a noticeable peak due to the J/^-neutron 
rescattering. The p — n rescattering peak is much larger than, and close enough to, 
the J j if>- neutron rescattering peak that it obscures the J/if> peak. Note also that the 
position of each of the peaks is simply given by the value of 6 n where the value of the 



corresponding |n_| is zero (see Fig. 4.35). 

The amplitudes F 2 b and F^, where the J /if) is produced on the neutron and 
then rescattering (of the neutron or J/if>, respectively) occurs on the proton, are 
much smaller than F 2a and F 3 b, and do not exhibit the well-defined peaks that F 2 b 
and F 3 b do. Fig. 4.39| shows the 8-fold electroproduction differential cross-section, 



Eq. (4.5.60), versus 6 n . Graphs are shown for 3 different values of the J/V'-neutron 
elastic scattering parameter Ay n : Ay n = 1.6 GeV~ 4 (which is the value determined 
in the experiment at SLAC), Ay n = 8.0 GeV~ A , and Ay n = 16 GeV~ A . It is seen that 
only if Ay n is of the order of 10 times as large as the previously measured value is there 
a noticeable peak due to the J /if}- neutron rescattering, for the "plus" kinematics. The 
p — n rescattering peak is much larger than, and close enough to, the J/i/'-neutron 
rescattering peak that it obscures the J /if) peak. For the "minus" kinematics, the 
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Figure 4.38: \F la + F 2a + F 3a \ 2 vs. 9 n , for 3 values of A Vn , for v = 10 GeF, t = 
—2 GeV^ 2 . The large peak on either side is due to p — n rescattering, while the small 
peak (or bump) on either side is due to J / -^-neutron rescattering. The dashed curve is 
only including the on-shell amplitude, while the solid curve includes both the on-shell 
and the off-shell amplitude. Here we show the "plus" and "minus" kinematics on the 
same graph by defining 6 n to be negative for the "minus" kinematics. 
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same statement holds; in addition, however, the size of the p — n rescattering peak 
varies (by ~ 40%) as the value of Ay n is varied. 

It's important to note that the peak due to the J/^-neutron rescattering isn't 



observable at lower energies. Fig. 4.40 shows the square of the total amplitude for 
photon energy of v = 9 GeV and t = —3, for Ay n = 10 x 1.6 GeV~ 4 . On this 
graph the peak due to p — n rescattering is visible, but there's no visible peak due to 
J I ?/>-neutron rescattering. 

4-6.2 Conclusion 

We have shown here the possibility of measuring the J/ -?/>-nucleon scattering amplitude 
for energies significantly smaller than the energy of the only existing data. If the 
total J/V'-nucleon cross-section a^ N at these energies is of the order of 2 — 3 times 
the previously measured value, then the differential cross-section as a function of 6 n 
should exhibit well-defined peaks corresponding to on-mass-shell p — n and J/ip — n 
rescattering, for virtual photon energy of v = 10 GeV and 4-momentum-transfer- 
squared t = (q — pv) 2 = —2 GeV 2 . However, at lower photon energy (below 9 GeV) 
the J/ip — n rescattering peak would not be distinguishable. As it is expected [33] 
that o J t lf N should increase as the energy decreases, it is not impossible that the lower 
energy cross-section could be larger than the measured value by a factor of ~ 2. 
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Figure 4.39: Electroproduction differential cross-section vs. 6 n , including all diagrams, 
for photon energy v = 10 GeV, Q 2 = 0.5 GeV 2 and t = -2 GeV 2 , for 3 values of 
A Vn : Solid curve: A y „ = 1.6 GeV- 4 . Dashed curve: A Vn = 8.0 GeV^ 4 . Dotted 
curve: Ay n = 16 GeV" 4 . In (a), the large peaks at 9 n ~ 45° and 9 n ~ —85° are due 
to p — n rescattering, while the small peaks (or bumps) at 9 n ~ 80° and 6>„ ~ —70° 
are due to J/^-neutron rescattering. (b) shows detail of left half of (a). 
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Figure 4.40: Amplitude squared vs. 9 n , including all diagrams, for Ay n = 16 GeV~ A , 
for photon energy v = 9 GeV and t = —3 GeV 2 . The peak or bump on either side 
is due to p — n rescattering, while the peak due to J/^-nucleon rescattering is not 
visible. 
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Appendix A 

GLAUBER THEORY RELATIONS FOR 2-BODY 

SCATTERING 

A.l Scattering amplitude and profile function 

Under the conditions for which the Glauber model is valid (high incident energy, small 
scattering angle, neglect of longitudinal momentum transfer and energy transfer) the 
elastic scattering amplitude /(q), where q = k — k' is the momentum transfer, can 
be written as the 2-dimensional Fourier transform of a function T(b) known as the 
"profile function": 

/(q) = — [ d 2 b e iqb r(b) (A.1.1) 



2vr „ 

where the vector b is perpendicular to the incident momentum k. If the scattering 
amplitude is known, then the profile function may be obtained by the inverse Fourier 
transform: 

r < b » = 2^/^^ b/ ( q )- < A ' L2 » 
The profile function T can be written in terms of a phase shift function xO 3 ) as 

r(b) = 1 - e ix{h) (A.1.3) 

and the scattering amplitude in terms of x is 

/(q) = ^ f d 2 b e iqb (l - e ix{h) ) (A.1.4) 



2tt . 

This shows the relation to the eikonal approximation in potential scattering: in the 



eikonal approximation, the scattering amplitude is given by Eq. (A.1.4) where the 
phase shift function is related to the potential V(r) by 

X (h) = --J V(b,z)dz. (A.1.5) 
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A. 2 Cross-sections 



The elastic, total, and inelastic cross-sections can be written in terms of the profile 
function. The differential cross-section is given by 

da 



dtt 



l/(q)l' 



(A.2.1) 



and so the cross-section for elastic scattering is 



a el = / dfi|/(q)| 5 



(A.2.2) 



For small angle scattering, the solid angle element dVt can be written as dVt = 
and so 

a * = hi ^ l/(q)|2 ' (A ' 2 - 3) 



which by using Eq. (A. 1.1) is 



0"eZ 



1 k 2 
k 2 (2tt) 2 



j d 2 q J d 2 bd 2 b' e iq - (b - b ' } r(b)r*(b') (A.2.4) 



which gives 



a el = / d 2 b\T(b)f 



(A.2.5) 



The optical theorem gives the total cross-section in terms of the forward scattering 
amplitude as 

An 4-tt h r 

a tot = y/m/(0) = — —Re J d 2 b T(b) (A.2.6) 



otot = 2 Re J d 2 b T(b) 



(A.2.7) 



Note that for /(0) pure imaginary (which holds to good approximation for high-energy 
hadron-hadron scattering) we have 



J d 2 bT(h) = l -a tot . 



(A.2.8) 
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Appendix B 



GLAUBER MODEL CROSS-SECTION FOR tt + A -> tt + A* 



The Glauber model calculation of the incoherent cross-section for h + A — > h + 
A* proceeds as follows [18]. Using the closure approximation (which neglects the 
dependence of the final states / of A on energy) the sum of the cross-sections for each 
final state is given by summing the squares of the amplitudes; then using closure on 
/ we obtain 



Assuming an independent particle model for the ground-state wavefunction 



<^|r to *(b)r^(bO|z>=i-[i-r(b)] A -[i-r*(bO] A +[i-r(b)-r*(bO+/3(b,bO] A 

(B.0.2) 

where f(b) = / d 3 r Pl (r)T(h - s) and /3(b,b') = / d 3 r Pl (r)T(h - s)T*(b' - s). And 



< i\v tot (b)\i >< i\rum >= i-[i-r(b)] A -[i-r*(b')] A +[i-r(b)] A [i-r*(b')] A . 

(B.0.3) 




i.e. |$i 



Yif=i Pi{ r j)i with J d 3 rpi(r) = 1 and p(r) = Api(r) we have 



also 
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Thus 



< i|r tot (b)r t * ot (b')|i > - < i\r tot (b)\i >< i\r* tot (b')\i > 

= [1 - r(b) - r (b') + /3(b, b')] A - [1 - r(b)] A [i - T* (b')] A 

= [(i - r(b))(i - r (b')) + p(b, b') - r(b)r* (b')] A - [(i - r(b))(i - r (b'))] A 

= [(i-r(b))(i-F(b'))] A 

+ - r(b))(i - F (b'))]- 4 - 1 (/?(b, b') - r(b)r*(b')) + •••-[(!- r(b))(i - r (b'))] A 
= a[(i- r(b))(i - r (b'^-^Cb.b') - r(b)r*(b')) + . . . 



The term proportional to r(b)r*(b') above only contributes significantly for q ~ 0. 
It's contribution to J2f^i l-^/^Q)! 2 i s 



(B.0.4) 



where the . . . represents higher powers of (/3(b,b ; ) — r(b)r*(b')). 




(B.0.5) 



The integral over b here is 




2tt 
ik 



/(q) 
(B.0.6) 



which gives 




(B.0.7) 



Without the factor e 2°" T ( S ) this would be proportional to the form factor of the 
nucleon density. Hence it is negligible for large q. 
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So neglecting that term, we are left with the term proportional to j3: 
(A)2 f #btftie*^\l-T(h)] A -\l-T{h')] A - x A/3(b,b') 

2,71 J 

= J d 2 sdzp(s, z){^-) 2 j dVb'e l ^ h - h '^e-^ T ^e-^ T ^T(b - s)T(b' - s) 

/^■<-» r(b - s) r( b '- s ) 



d 2 sdzp(s,z)e- aT{s \ — f 

2,71 

\f(q)\ 2 j d 2 sdzp(s,z)e~^ 



(B.0.8) 



where in the third line we've used the fact that T(b — s) and T*(b' — s) are sharply 
peaked at b = b' = s, to evaluate e~^ aT ^e~^ aT ^ h '^ at b = b' = s and pull it out of 
the integral over b. Thus we have for the single-scattering result for the incoherent 
cross-section: 

( rlrr\ C 

(B.0.9) 



da 
dVL 



incoh 



|/(q)| 2 | d 2 sdzp(s,z)e-^ 



The physical interpretation of this is that the projectile travels through the nucleus 
to the point (s, z) where it scatters elastically from a nucleon at that point, with 
momentum transfer q. It then continues along approximately in the z-direction out 
of the nucleus. The attenuation factor e~ aT ^ represents the probability that the 
projectile is not absorbed by the other nucleons, on its incident path and on its 
outgoing path. 

In [18] , Glauber derives an expression for the incoherent cross-section that includes 



all powers of (3. When this is expanded to first order in f3 it gives Eq. (B.0.9). 
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Appendix C 

SUMMING OVER ALL FINAL STATES OF RESIDUAL 

NUCLEUS 



In this appendix I derive the Glauber result, for tt + A — > ir+p+(A — 1)*, where we 
sum over all final states of the residual nucleus (not just one- hole states). The result 
is a multiple-scattering series, the terms representing one elastic scatter of the pion, 
two elastic scatters of the pion, etc., with the overall momentum transfer q shared 
between the individual scatterings. The first term (the single-scattering term) gives 



exactly the result obtained by summing over only one- hole final states, Eq. (2.4.26) 
which is as it should be. 

Starting from the general expression for the scattering amplitude for a particular 
final state / of the residual nucleus 

F fi = £ / d2b ^ h I • • • d 3r Ax;(n) 

2n J J (C.0.1) 

x 0A-i( r 2 ; • • • ,r A )0 j4 (r 1 ,r 2 , . . . ,r A )F tot (b, {r,}) 
and then squaring and summing over all / using the closure approximation gives: 

E l F /*! 2 = (^) 2 / d 2 be^ h d 2 b'e-^ W J dSdVidS d 3 r 3 . . . d 3 r AX ;(r 1 )x P (r' 1 ) 

(C.0.2) 

x 01(ri, r 2 , • • • , r A )(j) A (r l , r 2 , . . . , r A )r^ ot (b', r' 1; r 2 , . . . , r A )T tot (h, n, r 2 , . . . , r A ). 



(C.0.3) 



Now substituting 

.4 



l >* A (r f 1 ,T2,...,r A )<l> A (T 1 ,T2,...,T A ) = pii^Ti) JJpi^), (C.0.4) 

i=2 
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where in the shell model pi(r[,ri) = jr Yl,n=i ( t ) n{ r 'i) ( t ) n{ r i) with pi(r, r) = pi(r), and 
using orthogonality of the single-particle wavefunctions \ and 4> n , gives: 



k 



f 71 J J 



b'V 



(C.0.5) 



x 



/ d 3 r 2 d 3 r 3 . . . d z r APl {v 2 ) . . . Pl {r A ) f[(l - T bj - T* Vj + r 6 ,T 6 %.) (C.0.6) 
dVb'e 1 *^ J d 3 r 1 d 3 r' lX ;(r 1 )x P (r' 1 )p 1 (r' 1 ,r 1 )T bl T vv 



(C.0.7) 

x [1 -T(b) -T*(b') + f3(b,b')} A ~ 1 
where f(b) = / d 2 sdz Pl {s, z)r(b-s) and/3(b,b') = f d 2 sdzp 1 (s, z)r(b-s)r*(b'-s). 
This is to be compared to the expression obtained by squaring and only summing over 



one-hole states of the residual nucleus, which is Eq. (2.4.15) squared and summed: 

, k , 



•27T' 



d 2 bd 2 b'e^< h - h ^ / d'nd^MXpirdpM^bt r„ v g(b)g*(b') 



(C.0.8) 



where 



^(b)^(b') = [i - r(b)] A_1 [i - r*(b')] A_1 = [i - f (b) - r*(b') + f (b)r(b')] A -\ 

(C.0.9) 



Comparing the two expressions Eq. (C.0.7) and Eq. (C.0.8) we see that the difference 



lies in the last term inside the [. . .} A 1 . For the one-hole final states, it factorizes in 
b and b', while for the sum over all final states it does not. 



The result Eq. (C.0.7) can be expanded in powers of /3(b,b') to give a multiple 
scattering series: 

[1 - f (b) - f *(b') + /3(b, b')]^ 1 = [1 - T(b) - T*(b')] A - 1 

+ [1 - f (b) - f *(b')] A ~ 2 (A - l)/3(b, b') + . . . 

(C.0.10) 



(l + A/3(b,b') + ...) 



(C.0.11) 
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ff (b)0*(b')(i + 40(b,b') + ..O 



(C.0.12) 



using the approximation Eq. (2.4.19) (and in the large- A limit). The first term here 



thus contributes to J2j Wfi 

, k 4 



2 . 



'2vr' 



d 2 bd 2 b'e^ h - h ^ / C / 3 r 1 rfV lX ;(r 1 ) Xp (r' 1 )p(r' 1 ,r 1 )r fel r b , 1 ^(b)/(b / ) (C.0.13) 



which is exactly the same as the result obtained by only summing over one-hole states, 



Eq. (C.0.8 ). Keeping one power of (3 gives a term which corresponds to two scatterings 



of the pion, with the momentum transfer q shared between the two scatterings. 
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Appendix D 

LOW-ENERGY SCATTERING WAVEFUNCTIONS FOR 
J/^-NUCLEON SCATTERING 

The Schrodinger equation for the J/^-neutron system is: 



r 2 dr dr 2 



V le2 (r) = 2 f i,V(r)V k2 (r) (D.0.1) 



for I = (S-wave scattering). Normalization conventions: the wavefunction \I/ is 
related to u{r) and the phase shift S(k) by 

*( r ) = e «(fc)^fe) (D.0.2) 

This also defines our normalization of plane wave states: , 1 e lk r . 
The radial wavefunction u(r) satisfies 
d 2 u 

— + k 2 2 u = 2uV(r)u(r) = U(r)u(r) (D.0.3) 
dr 2 

We assume a finite range potential, so that V(r) is zero for r larger than some 



distance R. For r > R, where the potential is zero, the general solution to Eq. |D.0.3 
is 

u out (r) = — sin (k 2 r + 5{k 2 )) (D.0.4) 

which defines the S"-wave phase shift Sfa)- 

The phase-shift S(k) satisfies the following well-known properties [65] as k — » : 

1. for a repulsive potential, or an attractive potential that doesn't admit a bound 
state: 5 — > —ak as k — > 0; 



2. for an attractive potential which admits a single bound state: 5 — > n — ak as 
k ->■ 
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where a is a constant (the scattering length). Thus in the first case the zero-energy 
radial wavefunction for r > R is 

<ut(r) =r-a (D.0.5) 

and in the second case it is 

<*(r) = -(r - a), (D.0.6) 
while the zero-energy wavefunction ^„ ut is given in both cases by 

With these conventions, a > for either a repulsive potential or an attractive 
potential with a bound state, and a < for an attractive potential that doesn't 



admit a bound state (see Fig. 4.18). 



